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Abstract 

Fatigue  is  a  potential  failure  mode  of  submarine  pressure  hulls.  Australia's  new  Type  471 
submarines  will  be  constructed  from  a  new  high  strength  steel.  The  higher  design  loadings 
made  possible  by  the  higher  strength  make  fatigue  an  even  more  significant  consideration 
than  in  other  designs.  The  stress  analysis  of  submarine  pressure  hulls  under  hydrostatic 
loading  is  reviewed.  The  main  emphasis  is  on  analytical  methods  of  calculation  as  opposed 
to  numerical  (e.g.  finite  element)  methods,  and  on  the  calculation  of  the  fatigue  loading 
rather  than  the  calculation  of  buckling  pressure.  Because  of  the  emphasis  on  fatigue,  stress 
concentrations,  residual  stresses  and  the  interaction  between  residual  and  applied  stresses 
are  considered. 
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Analytical  Calculation  of 
Fatigue  Loading  of  Submarine  Hulls 


1.  Introduction 

The  following  introductory  remarks  arc  based  on  reviews  of  the  subject  by  Wenk  (1961), 
Kendrick  (1985)  and  Bushnell  (1985).  Submarine  pressure  hulls  are  the  world’s  largest 
pressure  vessels.  The  major  design  loading  is  the  external  hydrostatic  pressure  at  depth,  but  a 
submarine  hull  is  also  subject  to  hydrodynamic  stresses,  and  possibly  to  stresses  arising  from 
enemy  attack.  The  stress  analysis  of  the  hulls  is  a  well -developed  science  because  it  is 
important  to  produce  a  hull  of  minimum  weight  in  order  to  retain  sufficient  buoyancy  for 
propulsion  systems,  weapons  and  sensors.  Accurate  stress  analysis  results  in  precision 
design  and  the  confident  use  of  safety  factors  of  the  order  of  1 .5  or  less.  Accuracy  in  strength 
prediction  has  been  confirmed  by  tests  to  destruction  of  small  and  large  scale  models,  and  by 
strain  gauging  of  submarines. 

A  submarine  pressure  hull  typically  consists  of  a  cylindrical  centre  section  closed  at  the 
ends  with  spherical  or  torispherical  caps.  Some  examples  arc  shown  in  Figure  1 .  The  caps 
arc  usually  convex,  but  may  be  concave  (as  viewed  from  outside  the  submarine)  or  even  flat 
as  in  many  Kockums  designs.  The  cross  section  may  be  tapered  towards  the  ends  by  the 
insertion  of  conical  sections.  More  general  shapes  are  also  possible,  including  figure-of-eight 
or  even  triple-hulled  cross  sections,  or  more  fish-shaped  single  hull  configurations.  A  recent 
Italian  design  consists  of  numerous  toroids  welded  together  (Compton-Hall,  1989). 

The  challenge  for  stress  analysts  is  to  calculate  the  extent  of  stress  concentration  at 
transitions  (c.g.  between  the  end  caps  and  remainder  of  the  hull,  cone-cylinder  transitions),  at 
hull  penetrations,  and  at  other  structural  discontinuities  such  as  bulkhead  locations. 


(a)  German  Type  XXI  U-boat 


(g)  Dutch  Zwaardvis  class 


|C«y»>fnootc*«tO't  M*»t 


(h)  British  Type  2400 


(i)  Italian  3GST9  midget  submarine 
Figure  1  (continued) 
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2.  Current  Relevance 


This  review  is  concerned  with  the  stress  analysis  of  Australia’s  new  Type  471  submarines  for 
the  purpose  of  assessing  susceptibility  to  fatigue  and  corrosion  fatigue.  The  pressure  hull  of 
the  Type  471  is  being  constructed  from  a  new,  high  strength,  qucnched-and-tempered  steel. 
The  susceptibility  of  steels  to  corrosion  fatigue  generally  increases  as  their  strength  increases, 
only  partly  because  the  higher  strength  allows  design  lor  higher  stresses  (de  Hart,  1969; 
Thomson  and  Christopher,  1973).  In  addition,  a  submarine  can  expect  a  fatigue  loading  of 
the  order  of  1()4  cycles  during  its  lifetime  due  to  diving  to  depth.  This  puts  it  in  the  regime  of 
low  cycle  or  medium  cycle  fatigue,  where  increased  yield  strength  may  decrease  the  S-N 
fatigue  life  (Stout  and  Pensc,  1965).  Therefore  the  Type  471  is  potentially  much  more 
susceptible  to  fatigue  and  corrosion  fatigue  than  the  Oberon 

In  this  context,  stress  analysis  is  not  required  for  calculation  of  collapse  depth  or 
buckling  mode.  With  the  Type  471  design  akcady  established,  the  requirement  is  for 
calculation  of  stresses  in  the  hull  as  a  function  of  depth  down  to  about  two  thirds  of  collapse 
depth  (assuming  a  safely  factor  of  1 .5).  W  thin  this  range  an  elastic  stress  analysis  is 
sufficient  to  obtain  the  necessary  data  for  fatigue  and  corrosion  fatigue  analysis.  The 
combination  of  the  expected  mission  profile  and  service  life  of  the  Type  471  submarine 
together  with  the  elastic  stress  analysis,  allows  a  fatigue  loading  spectrum  to  be  predicted. 

In  tin's  report  the  emphasis  is  on  analytical  solutions  involving  hand  or  simple  computer 
calculations.  This  provides  a  qualitative  picture  of  the  overall  behaviour  and  provides 
calibration  points  for  more  sophisticated  computer  stress  analysis.  Axisymmetric  stress 
analysis  has  become  a  routine  process  using  codes  such  as  ASSSAI  and  BOSOR4  (Roberts 
and  Smith,  1 988). 


3.  Simple  Structures 


In  this  section  some  simple  structures,  idealisations  of  a  submarine  hull,  arc  examined.  In 
subsequent  sections  more  realistic  configurations  are  considered.  Consider  first  an  infinitely 
long,  circular  lube  subjected  to  uniform  external  pressure.  The  circumferential  (hoop)  stress 
is  (Young,  1989): 

a()  =  -  pR  /h  (1) 

and  the  reduction  in  radius  is: 

AR  =  pR:/Eh  (2) 


where  p  is  the  pressure,  R  the  radius  and  h  the  wall  thickness  of  the  cylinder. 

A  finite-length  cylinder  experiences  similar  hoop  stress  but  is  also  subject  to  a 
longitudinal  stress  of: 


M 


a/  =  ct0/2  =  -  pR  /  2h 


0) 


due  to  pressure  on  the  end  closures.  The  basic  stress  slate  of  a  submarine  hull  is  there  lore 
one  of  biaxial  compression  with  the  longitudinal  stress  half  that  in  the  circumferential 
direction.  Stress  concentrations  in  a  real  hull  can  modify  this  significantly  in  local  regions 
and  even  result  in  regions  of  tensile  stress. 

The  reduction  in  radius  for  the  finite  length  cylinder  is: 


AR  =  pR2  (1  -  v/2)/Eh  (4) 

Minimum-weight-design  of  a  hull  necessitates  the  use  of  a  stiffened  cylinder  rather  than  a 
simple  shell.  Because  of  the  large  circumferential  stress  component,  the  stiffeners  are  in  the 
form  of  circumferential  rings  welded  to  cither  the  inside  or  the  outside  ol  the  hull  plating. 

The  stiffening  rings  (or  "frames")  usually  have  a  T  cross  section.  Under  pressure  the  shell 
tends  to  bend  inward  between  the  stiffeners  forming  a  concertina  shape.  Bending  stresses  are 
introduced,  and  the  stresses  are  no  longer  unifonn  through  the  thickness  of  the  plating  as  they 
are  in  the  case  of  a  simple  tube.  The  schematic  geometry  is  shown  in  Figure  2. 


4.  Yielding ,  Buckling  and  Submarine  Design 

Failure  of  a  pressure  hull  at  great  depth  may  be  due  to  yielding  of  the  hull  steel,  or  due  to 

some  form  of  buckling.  The  design  determines  which  mode  of  failure  actually  lakes  place. 

The  four  primary  modes  of  failure  to  be  considered  (Wenk,  1961;  Kendrick,  1985)  are: 

(i)  Yielding  of  the  shell  between  the  rings,  characterised  by  an  axi-symmclric 
concertina-shaped  pleat.  This  is  also  the  elastic  deformation  mode  prior  to  collapse. 
Yielding  and  other  failure  modes  arc  illustrated  in  Figure  3. 

(ii)  Inter-stiffener  buckling  ("local  instability",  "shell  buckling")  in  which  lobes  develop 
between  the  rings.  The  rings  are  at  nodes  of  the  buckling  mode,  and  there  is  usually  a 
large  number  of  lobes. 

(iii)  Stiffener  tripping  ("local  stiffener  instability")  in  which  a  ring  fails  by  buckling. 

(iv)  General  instability  ("overall  collapse")  in  which  rings  and  shell  lail  simultaneously 
from  one  end  of  the  stiffened  cylinder  to  the  other,  with  two  to  live  circumferential 
lobes.  In  practice,  general  instability  is  probably  precipitated  by  initial  yielding  or 
stiffener  tripping. 
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Figure  2:  Schematic  diagram  showing  deflection  of  a  ring-stiffened  cylinder  under 
hydrostatic  pressure  to  form  a  concertina  shape. 


Lobes  Pleat 


Figure  3:  Modes  of  collapse  of  a  ring-stiffened  cylinder  (from  Roberts  and  Smith,  1988). 


Pressure  hull  failure  modes  mirror  the  failure  modes  of  a  column  under  compression 
loading.  Short,  thick  columns  fail  by  yielding,  as  do  hulls  with  thick  shells  and  closely 
spaced  stiffeners.  Slender  columns  and  thin  shells  with  widely  spaced  stiffeners  fail  by 
buckling.  I  he  dimensions  of  the  stiffeners  determine  whether  shell  buckling  is  local  (heavy 
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stiffeners)  or  overall  (light  stiffeners),  or  whether  tripping  occurs  (very  light  stiffeners  or 
those  with  thin  webs). 

To  illustrate  the  general  principles  of  buckling,  consider  a  simple  column.  A 
slender  column  under  axial  compression  remains  straight  up  to  a  critical  load,  Pcr,  where 
buckling  may  occur.  Below  Pcr  the  column  is  in  stable  equilibrium.  Above  Ptr  two 
equilibrium  positions  arc  possible  (one  straight,  the  other  buckled).  Even  a  small 
perturbation  can  cause  the  transition  from  straight  to  buckled  configuration.  This  is  known 
as  bifurcation  buckling. 

For  a  pin-ended  column  die  critical  buckling  load  is  (Popov,  1978): 

P.,  =  "2e'x'LJ  (5) 

where  E  is  Young’s  modulus,  I  is  the  moment  of  inertia  and  L  the  length  of  the  column. 
The  corresponding  critical  compressive  stress  in  the  column  is: 

CT  =  P  /A  =  tc2E/X2  (6) 

cr  cr  v 


where  A  is  cross  section  and  X  the  slenderness  ratio,  \I  /  L2A. 
Above  the  critical  load  or  stress,  the  buckled  mode 

w  =  w()  sin  (k/./L) 


(7) 


becomes  possible  in  addition  to  the  straight  column  (axisymmetric)  mode  w  =  0,  where  w  is 
the  sideways  deflection  and  z  the  coordinate  in  the  axial  direction. 

If  ocr  exceeds  the  yield  stress  of  the  material,  oy,  then  yielding  occurs  rather  than 
buckling.  As  yielding  occurs  and  work  hardening  begins,  the  effective  stillness  is  the  tangent 
modulus  E(  rather  than  the  elastic  modulus  E.  Therefore  bifurcation  buckling  may  occur  at 
the  point  of  yielding  or  after  some  plastic  deformation.  These  ideas  are  summarised  in 
Figure  4  (from  Wenk,  1961).  Parameters  analogous  to  ctcr  and  X  for  stiffened  cylindrical 
shells  arc  identified  and  test  results  compared  with  equations  analogous  to  equation  (6). 

True  bifurcation  buckling  occurs  only  in  unrealistic,  ideal  situations  where  the  buckling 
deformation  mode  is  orthogonal  to  the  clastic  deflection  mode  (i.c.  buckling  mode  amplitude 
equals  zero  regardless  of  the  extent  of  elastic  deflection;  Bushncll,  1985).  In  practice, 
geometrical  imperfections  and/or  loading  eccentricity  arc  present,  so  that  the  buckling  mode 
develops  continuously  from  the  initial  irregularities  as  load  is  increased.  Non-linear, 
clastic-plastic  collapse  (under  rising  load)  or  geometrical  "snap-through"  (under  falling  load) 
arc  the  only  buckling  modes  observed.  Increasing  eccentricity  or  initial  imperfection 
amplitude  decreases  the  buckling  load  as  shown  in  Figure  5. 
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Yielding 

j  Plastic  Instability 
I  Elastic  Instabilit 


Slenderness  X  =  V2RLf/n /cj/E 


Figure  4:  Experimental  data  points  showing  the  effect  of  slenderness  ratio  X  on  the  failure 
mode  and  collapse  pressure  pc  for  ring-stiffened  cylinders.  Elastic  buckling  occurs  for 
X>  l .2  (from  Wenk,  1961). 


IS 


(a)  Radial  deflection  of  an  imperfect  cylinder  on  increasing  pressure  from  0  to  SU  psi. 


Ordinary  Buckling  Snap  Through 


Displacement  Displacement 


(b)  Schematic  load-displacement  curves  for  ordinary  buckling  ( column )  and  snap  through 
(spherical  shell). 

Figure  5:  Effect  of  geometrical  imperfections  (or  eccentricity  e)  on  the  buckling  of  a 
cylinder,  column  and  sphere  (from  Wenk,  1961). 


16 


The  calculation  of  the  critical  pressure  (or  collapse  depth)  is  a  primary  concern  during 
the  design  of  a  submarine  hull.  The  usual  approach  is  to  consider  each  possible  failure  mode 
separately.  It  is  possible  to  predict  pcr  more  accurately  for  some  modes  than  for  others,  so 
different  safety  factors  arc  applied  for  the  different  modes  (for  example,  overall  collapse  is 
more  sensitive  to  imperfections  than  other  modes).  The  usual  procedure  (Kendrick,  1 985)  is 
to  proportion  the  shell  plating  thickness  and  stiffener  dimensions  so  that  collapse  at  pressures 
in  excess  of  the  design  pressure  would  occur  by  inter-stiffener  collapse  (with  a  safety  factor  of 
1.5).  Premature  failure  by  overall  collapse  or  local  stiffener  instability  are  avoided  by 
applying  safety  factors  of  2  and  4  respectively.  Part  of  the  reason  for  applying  different 
safety  factors  is  to  avoid  interaction  between  different  failure  modes. 

In  this  report  buckling  is  not  a  major  concern.  Rather  the  interest  is  in  stress  analysis  at 
external  pressures  encountered  in  normal  service,  well  below  the  critical  buckling  pressure. 
However,  considerations  of  buckling  largely  determine  the  submarine  design,  and  therefore 
indirectly  affect  the  stresses  and  fatigue  loading  encountered  in  normal  service. 


5.  Elastic  Stress  Analysis  of  Submarine  Hulls 

5.1  Introduction 

Much  of  the  literature  on  the  stress  analysis  of  ring-stiffened  cylinders,  and  other  idealisations 
of  submarine  pressure  hulls,  is  concerned  with  the  prediction  of  buckling  pressure.  The 
calculation  of  clastic  stresses  in  the  hull  is  a  much  simpler  task. 

The  early  work  was  carried  out  by  von  Sanden  and  Gunther  (1920)  and  von  Mises  (1929) 
in  Germany  at  the  time  of  World  War  I.  This  was  extended  after  the  war  by  various 
researchers  at  the  David  Taylor  Model  Basin,  as  reviewed  by  Wenk  (1961).  More  recent 
reviews  arc  those  of  Kendrick  (1970,  1985)  and  Bushncll  (1985). 

Work  prior  to  about  1955  emphasises  analytical  solution  procedures  and  is  restricted  to 
fairly  simple  (but  not  trivial)  geometries.  A  wide  variety  of  general  and  more  specific 
computer  programs  is  now  available  for  the  stress  analysis  of  structures  including  pressure 
vessels.  The  codes  dealing  with  axisymmctric  structures  are  most  applicable  to  the  stress 
analysis  of  submarine  hulls. 


5.2  Illustrative  Example 

In  order  to  gain  a  qualitative  picture  of  hull  stresses,  consider  the  "typical"  structure  shown  in 
Figure  6.  The  stress  distributions  at  various  locations  obtained  using  NCRE’s  PEGASUS 
computer  code  (Kendrick,  1964)  arc  shown  in  Figure  7.  Away  from  stress  concentrations 
other  than  uniformly  spaced  stiffeners,  the  general  features  of  the  stress  distribution  in  the 
cylindrical  sections  arc: 
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28  AT  45 


Figure  6:  Example  of  a  typical  submarine  pressure  hull  ( Kendrick ,  1 964 ). 


(a)  Stresses  near  dome-cylinder  junction 

Figure  7:  Stresses  in  the  hull  plating  for  the  structure  shown  in  Figure  6. 
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(b)  Stresses  near  cone-cylinder  junction. 


Figure  7  (continued) 


(i)  The  stress  levels  vary  periodically  on  passing  from  frame  to  frame. 

(ii)  The  average  longitudinal  stress  (i.e.  average  of  stresses  at  interior  and  exterior  surfaces) 
is  approximately  oy  =  -pR/2h  and  the  average  circumferential  stress  is  slightly  less  than 
Oq  =  -pR/h  due  to  the  stiffening  effect  of  the  rings  (oz  and  Oq  arc  the  stresses  for  a 
uniform  cylinder  given  in  Section  3.  For  a  cylinder  of  radius  3  m  and  thickness 

25  mm,  oy  =  -60  MPa  and  Oq  =  -120  MPa  at  100  m  depth). 

(iii)  Stresses  vary  in  the  longitudinal  direction  due  to  shell  bending  over  the  rings.  The 
longitudinal  stresses  show  the  most  variation,  with  the  external  and  internal  surface 
stresses  mirror  images  of  each  other  (reflected  in  the  mean  stress  -pR/2h). 

(iv)  The  circumferential  stresses  arc  similar  at  both  surfaces,  but  mimic  the  variation  of  the 
corresponding  longitudinal  stress  because  of  the  Poisson  effect. 

Bulkheads  and  section  transitions  may  cause  large  stress  concentrations,  although  their 
influence  docs  not  spread  very  far  into  the  adjacent  structure.  The  stress  concentrations 
occur  because  of  differences  in  the  radial  stiffness  of  adjacent  sections,  and  the  resultant 
bending  of  the  plating  to  accommodate  the  different  deflections  under  external  pressure.  The 
stress  concentrations  at: 

(i)  A  dome-cylinder  junction 

(ii)  A  cone-cylinder  junction 

(iii)  A  bulkhead  in  a  cylindrical  section 

arc  shown  in  Figure  7(a),  (b)  and  (c)  respectively.  The  stress  concentration  at  the 

dome-cylinder  junction  is  very  small  compared  with  the  other  two. 

The  stress  concentration  factors  can  be  reduced  by  curved  transitions  and  other  methods 
of  matching  the  radial  stiffness  across  the  structural  discontinuity.  This  may  be  done  by 
locally  modifying  the  shell  thickness,  frame  depth  or  frame  spacing. 


5.3  Symbols 

The  symbols  listed  below  are  used  in  the  calculation  of  the  hull  plating  stresses  in  following 
sections  of  this  report.  Some  arc  illustrated  in  Figure  8. 

0  circumferential  coordinate 

z  axial  coordinate  (measured  from  mid-way  between  frames) 

r  radial  coordinate  (measured  from  cylinder  axis) 

p  radial  coordinate  (measured  from  shell  mid-thickness) 
u  axial  displacement  of  mid-surface  of  shell 

v  circumferential  displacement 

w  radial  displacement  (positive  outwards) 

R  shell  radius  at  mid  thickness 

Rf  radius  of  centroid  of  ring  stiffener 


h  shell  thickness 

Lf  ring  stiffener  spacing 

b  ring  faying  flange  width 

L  unsupported  length  of  shell  (L  =  Lf  -  b) 

Af  cross  sectional  area  of  ring  stiffener 
p  hydrostatic  pressure  (positive  if  external) 

E  Young’s  modulus  of  shell  and  rings 
v  Poisson’s  ratio  of  shell  and  rings 
o  stress 

i  strain 

M  bending  moment 

Q  radial  shear  force 

N  membrane  force  (stress  resultants) 

I  h3  /  [  12  (1  -  v2)] 

D  flexural  rigidity  of  shell  (D  =  El) 

I,  moment  of  inertia  of  ring  stiffener 
a  fourth  root  of  3(  1  -  v2)  /  (hR)2 

K  frame  rigidity  constant  (total  load  on  ring  per  unit  radial  deflection  per 
unit  circumferential  length) 


Figure  8;  Symbols  used  to  describe  the  geometry  of  uniformly-framed,  ring-stiffened 
cylinders. 
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5.4  Uniformly  Framed  Cylinders 


The  first  complete  stress  analysis  of  a  perfectly  circular,  transversely-framed  cylinder  was 
carried  out  by  von  Sandcn  and  Gunther  (1920,  1921).  It  is  the  basis  of  subsequent  analytical 
solutions  of  the  problem,  both  in  the  formulation  of  the  differential  equation  and  in  the 
method  of  solution.  They  found  that  stresses  in  the  shell  plating  are  strongly  influenced  by 
ring-stiffener  cross  section  and  spacing. 

The  mathematical  analysis  considers  an  element  of  the  shell  consisting  of  a  longitudinal 
strip  from  one  stiffener  to  the  next.  The  strip  is  analysed  as  a  beam  supported  by  contiguous 
shell  elements,  with  its  ends  resting  on  an  elastic  foundation  (the  stiffening  rings).  Therefore 
the  governing  differential  equation  is  similar  to  that  for  pure  bending  of  a  beam: 

EIx  (d2w/dz2)  =  M  (8) 


The  full  derivation  of  the  differential  equation  is  given  by  Timoshenko  and 
Woinowsky-Kricgcr  (1959),  for  example. 

A  preliminary  analysis  was  carried  out  ignoring  the  effect  of  pressure  on  the  end  closures 
of  the  cylinder,  allowing  free  movement  in  the  axial  direction  so  o?  =  0.  The  resulting 
boundary  value  problem  and  its  solution  arc  presented  in  Appendix  1 .  Two  of  the  boundary 
conditions  are  obtained  from  symmetry  considerations,  and  the  third  is  obtained  by 
considering  the  radial  force  balance  at  the  stiffeners  (external  pressure  plus  shell  shear  force 
equals  ring  reaction  force). 

Von  Sanden  and  Gunther  also  modified  their  solution  to  partially  take  end  pressure  into 
account.  The  changes  arc  listed  in  Appendix  2.  The  longitudinal  membrane  stress  becomes 
-pR/2h  but  other  stresses  are  unchanged. 

Viterbo  (1930)  pointed  out  that  if  Af  is  set  to  zero  (i.e.  a  cylinder  without  ring 
stiffeners),  the  von  Sanden-Gunther  solution  predicts  a  shell  load  of  less  than  p  per  unit  area 
(the  applied  pressure).  The  inconsistency  is  removed  by  dividing  the  (radial)  pressure  into 
components  pv/2  and  p  -  pv/2,  and  noting  that  the  pv/2  component  together  with  the  axial 
load  due  to  end  pressure  produces  (approximately)  no  deformation  of  the  stiffening  rings  (sec 
Young,  1989).  This  leads  to  a  modification  of  the  third  boundary  condition  and  a  new 
solution  as  shown  in  Appendix  3. 

Salerno  and  Pulos  (1951)  further  refined  the  solution  by  incorporating  the  full  effect  of 
the  end  pressure.  The  Viterbo  boundary  conditions  still  apply  but  the  differential  equation 
includes  a  non-linear,  beam-column  term  (Appendix  4).  Solution  of  this  problem  is  along  the 
same  lines  as  the  von  Sanden-Gunther  solution  (sec  also  Moshaiov  and  Joclson,  1989)  but  it  is 
more  tedious  and  results  in  more  complicated  expressions  for  shell  deflection  and  stresses. 
However  it  is  regarded  as  the  exact  solution  for  a  uniform,  ring-stiffened  cylinder  under 
external  hydrostatic  pressure.  Deflection  and  stresses  are  no  longer  linearly  related  to  the 
pressure. 

The  process  of  solution  of  the  complete  Salcmo-Pulos  problem  has  been  simplified  for 
stresses  at  mid-bay  (at  z  =  0,  midway  between  the  stiffeners)  by  the  graphical  solutions  of 
Krcnzkc  and  Short  (1959).  Wilson  (1956  a,  b,  c)  used  two  methods  to  obtain  accurate 
approximations  to  the  Salcmo-Pulos  solutions.  One  method  assumes  that  tnc  shell  deflection 
as  a  function  of  pressure  takes  the  form  of  a  power  scries  in  p.  The  other  assumes  a  Fourier 
scries  form. 


22 


The  power  series  (asymptotic)  solution  is  described  in  Appendix  5.  For  various  test 
cases  studied  by  Wilson,  the  asymptotic  solution  is  within  5%  of  the  Salemo-Pulos  solution. 
It  is  regarded  as  satisfactory  provided  the  linear  term  pw,  is  much  larger  than  the  correction 
term  p2w?  (where  the  total  deflection  is  w  =  pwj+p2w2). 

The  Fourier  scries  method  (Appendix  6)  allows  another  small  refinement  to  the 
governing  differential  equation.  In  practice,  only  the  first  few  terms  in  the  Fourier  series  arc 
required  to  obtain  an  accurate  approximation  to  the  full  Salemo-Pulos  solution.  Wilson  also 
used  slightly  modified  expressions  for  the  frame  rigidity  constant,  K,  and  the  shell  stresses. 
The  stress  equations  are  based  on  simple  elasticity  equations  for  plane  stress  such  as: 

axx  =  E(e  +V€  )/(l-v2)  (9) 

xx  xx  yy 

Wilson  also  considered  the  effects  of  factors  such  as  out-of-circularity  on  shell  stresses. 
These  effects  arc  considered  in  Section  5.10. 


5.5  Computer  Program  for  Uniformly  Framed  Cylinders 

A  computer  program  has  been  written  at  MRL  to  implement  the: 

(i)  Von  Sandcn  and  Gunther 

(ii)  Viterbo 

(iii)  Wilson  power  series 

(iv)  Wilson  Fourier  series 

solutions  for  stresses  in  the  shell  of  a  uniformly  framed  cylinder  subject  to  external  pressure. 
Shell  displacement  and  curvature,  and  longitudinal  and  circumferential  stresses  in  the  shell  at 
mid-thickness  and  both  surfaces  arc  calculated  (shell  buckling  is  not  considered). 


5.6  Stress  Concentrations  due  to  Axisymmetric  Structure 

To  this  point  only  uniformly-framed  cylinders  have  been  analysed.  In  this  section, 
consideration  is  given  to  some  departures  from  uniform  framing  such  as  bulkheads,  heavy 
stiffeners,  transitions  (c.g.  cylinder-cone  junctions)  and  junction  of  a  cylinder  with  end 
closures. 

The  analysis  is  based  on  that  presented  by  Gill  (1970)  for  infinite  cylinders  subject  to 
edge  loading.  For  asymmetric  loading  of  a  shell  element  of  length  dz  and  width  r  d0  the 
stress  resultants  arc  as  shown  in  Figure  9  and  the  equilibrium  equations  arc  as  follows: 
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Figure  9:  Stress  resultants  for  asymmetric  loading  of  an  element  dzxr  dQ  of  a  cylindrical 
shell. 
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If  the  loading  is  axially  symmetric  (c.g.  hydrostatic  pressure),  pg  =  0  and 
N,e  =  Ngz  =  M/0  =  Qg  =  0  so  the  equations  reduce  to: 


(dN^/dz)  +  pz .  =  0 
r  (dQz/dz)  +  N0  -  rpr  =  0 
(dM/lz)  -  Qz .  =  0 


(11) 


For  linear  clastic  material,  and  considering  only 
z  =  const.),  die  following  simplifications  apply: 

Pr  =  P,  =  0 

N/  =  0 
N0  =  Ehw/r 

M;  =  D(d2w/dz2) 

Substituting  into  the  equilibrium  equations  (11) 

El  (d4w  /  dz.4)  +  Ehw  /  R2  =  0 

This  is  the  differential  equation  for  a  circular  cylinder  subjected  to  rotationally  symmetric 
edge  loading.  It  is  equivalent  to  the  initial  Van  Sandcn  and  Gunther  analysis  with  p  =  0 
(Appendix  1).  The  general  solution  is  as  in  Appendix  1  with  p  =  0,  or  equivalently,  and  more 
conveniently  in  this  case,  is: 

w  =  exp  (cl/.)  |kj  cos  (az)  +  k.,  sin  (az)] 

(14) 

+  exp  (-az)  (k3  cos  (az)  +  k4  sin  (az)] 


shear  (Qz)  and  moment  (Mz)  edge  loading  (at 


(12) 


gives: 


(13) 


Consider  a  semi-infinite  cylinder  subjected  to  symmetric  edge  loadings  Mo  and  Qo  at  z  =  0 
(Figure  10).  Because  the  effect  of  the  edge  loading  must  decay  with  increasing  z,  the 
coefficients  k,  and  k.,  arc  zero.  The  coefficients  k3  and  k4  arc  determined  by  the  boundary 
conditions  at  the  free  edge  (z  =  0). 

Due  to  the  exponential  term,  the  radial  deflection  w  is  significant  only  for  az  <  rc.  Fora 
5  m  diameter  steel  pressure  vessel  with  walls  25  mm  thick,  this  condition  is  equivak.it  to  z  < 
0.71  m  (this  distance  is  of  the  order  of'lhc  spacing  of  submarine  hull  ring  stiffeners).  Thus  the 
effects  of  the  edge  loading  arc  significant  only  in  one  bay. 

This  is  an  important  simplification  in  the  analysis  of  hull  stresses.  A  stress 
concentration  such  as  a  bulkhead  or  heavy  stiffener  imposes  radially  symmetric  loads  on  the 
edge  of  die  adjacent  hull  cylinder.  The  resulting  deflections  will  usually  only  be  significant 
in  the  hull  plating  between  the  stress  concentration  and  the  nearest  stiffening  ring.  However, 
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if  two  stress  concentration  points  are  located  closer  together  than  nta,  then  the  exp(az)  term 
in  equation  (14)  cannot  be  deleted  and  the  interference  between  the  two  must  be  considered. 


Figure  10:  Edge  loadings  Qo  and  Mo  on  a  cylindrical  shell  produce  deflection  wo  and  slope 
tan  (9 J  (from  Gill.  1970). 


Table  1  (Gill,  1970)  shows  solutions  of  equation  (14)  for  various  boundary  conditions  (edge 
loadings).  In  the  case  of  a  submarine  hull,  the  deflections  due  to  hydrostatic  pressure  must  be 
added  to  obtain  the  complete  solution.  This  process  is  described  for  various  types  of  stress 
concentrator  in  the  following  sections. 

5.6.1  Cylinder -Hemisphere  Junction 

A  well-studied  situation  giving  rise  to  a  stress  concentration  is  the  closure  of  a  cylinder  with  a 
hemispherical  cap.  In  the  case  of  equal  shell  thicknesses,  the  stresses  and  radial 
displacements  of  simple  cylindrical  and  spherical  shells  due  to  external  pressure  are: 
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Table  1:  Edge  Solutions  for  a  Long  Cylindrical  Shell 
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Cylinder: 

Oq  =  -pR/h 

0) 

at  =  -  pR  /  2h 

(3) 

w,  =  -  pR2  (2  -  v)  /  2Eh 

(4) 

Sphere: 

CTe  =  °v  =  -PR/2h 

(15) 

8,  =  -  pR2  (1  -  v)  /  2Eh 

(16) 

The  deflections  are  not  equal,  so  equal  and  opposite  shear  forces  and  bending  moments  must 
exist  at  a  junction  between  cylinder  and  hemisphere  in  order  to  maintain  continuity.  The 
edge  loading  solutions  for  a  cylinder  shown  in  Table  1  and  similar  results  for  a  sphere  shown 
in  Table  2  allow  the  shear  force  Qo  and  bending  moment  Mo  to  be  obtained  by  equating  the 
displacements  wo  and  5o  and  slopes  dw/dz  and  tanjfl^  on  either  side  of  the  junction.  The 
procedure  for  calculating  the  total  stress  in  the  cylindrical  and  hemispherical  shell  involves 
superposition  of  the  membrane  and  edge  bending  solutions,  and  is  known  as  "discontinuity 
analysis"  (Cloud,  1970). 

Gill  (1970)  gives  details  of  the  calculations.  Referring  to  Figure  1 1  and  Tables  1  and  2 
(noting  p  =  tt/2  in  this  case),  the  edge  displacements  and  slopes  due  to  edge  forces  Qo  and  Mo 
are: 


Cylinder: 

wo  =  2Q0  aR2  /  Eh  +  2M0  a2R2  /  Eh 

(17) 

9o  =  -2Qoa2R2/Eh-4Moa3R2/Eh 

(18) 

Sphere: 

5o  =  2Mq a2R2 /  Eh - 2Qo aR2  /Eh 

(19) 

eo  =  4Mo  “3r2  /  Eh  ’  2Q0  a2R2 1  Eh 

(20) 

Matching  the  deflections  and  slopes  across  the  discontinuity  due  to  the  combination  of 
hydrostatic  pressure  with  the  edge  loadings  therefore  requires: 

w1  +  wo  =  51  +  60  (21) 

9o  =  -2Qo  a2R2/Eh  -  4Mo  a3R2/Eh  =  4Mo  ct3R2/Eh  -  2Qo  a2R2/Eh 

(22) 


These  allow  M  and  Q  to  be  obtained: 

O 

M  =  0 

O 

Qo  =  p/8a 


(23) 

(24) 
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Table  2  Edge  Bending  Solutions  for  a  Spherical  Shell  (4)  >  30°) 


Geometry 


Stress  Resultant 

A  exp(-aRy)  sm(«Rvy  +  y) 

N'0  -  Q0  col(P  -  v) 

N0  -  V2aR  A  cxp(-aRv|/)  sin(aRy+y-n/1) 

+(I/N2a)  A  exp(  aRH')  sin(aRv+y+ic/4) 

M0  vM0 

A  exp  (  aRy)  sin(aR\(/  +  y) 

-  Qty  coi(P  +  y) 

+  V2aR  A  exp(-aR\)/)  sin(aRv*y-rt/4) 

-(1  /V2a)  A  exp(-aRiff)  sin(aRif+y+n/4) 

vM. 

<1> 

Rotation 

9 

(2a2R2/Eh)  A  cxp(-aRy)  cos(aRiy+y) 

(2a2R2/Eh)  A  exp(-aRv)  cos(aR\|/+y) 

Constants  for  Boundary  Conditions* 

Mq  only 

A  =  2aM  y=0 

o  * 

A  =  -2aM  v=0 

o  1 

Qo°nly 

A  =  \2  Qo  sin  P  y  =  -n/4 

A  =  -V2  Qq  sin  P  y  =  -rc/4 

Edge  Rotation 

0o 

(4a3R2/Eh)MQ  +  (2a2R2/Eh)QQ  sin  P 

-4(a3R2/Eh)M0  -  (2a2R2/Eh)Q0  sin  P 

Edge  Deflection 

8o 

(2a2R2/Eh)MQ  sin  P  +  (2aR2/Eh)Q0  sin2  P 

(2a2R2/Eh)MQ  sin  P  +  (2aR2/Eh)Qo  sin2  P 

The  values  of  A  and  y  are  found  for  M0  and  Q0  applied  separately  and  added  before 
substituting  into  the  expressions  for  the  stress  resultants  and  rotation. 
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Figure  1  / :  Deflections  and  slopes  at  a  hemisphere-cylinder  junction  (a)  due  to  pressure  and 
(b)  due  to  edge  loadings  Qo  and  Mo. 


The  total  deflections,  due  to  Q0  (obtained  trom  Tables  1  and  2)  combined  with  those  due  to 
hydrostatic  pressure  are: 

w  +  Wj  =  (pR2/4Eh)  exp(-az)  cos(az)  -  pR2  (2-v)  /  2Eh  (25) 

8  +  8,  =  -(pR2/4Eh)  exp(-Ra\)/)  cos(Rav|/)  cos(\)/)  -  pR2  (1-v)  /  2Eh  (26) 

Knowing  the  deflections,  the  stresses  can  now  be  calculated  as  in  Appendices  1  to  6.  Results 
are  shown  in  Figure  12  (Gill,  1970).  Note  that  the  bending  of  the  sphere  and  cylinder  near 
the  junction  increases  the  longitudinal/meridional  stress  from  -pR/2h  to  about  -1.3  pR/2h  at 
one  surface  (and  reduces  it  to  -0.7  pR/2h  at  the  other).  The  hemisphere  is  stiffer  than  the 
cylinder  in  the  radial  direction,  so  the  radial  deflection  of  the  cylinder  is  reduced  at  the 
junction.  This  leads  to  reduced  circumferential  stress  at  that  point  for  the  cylinder  (and 
increased  stress  for  the  sphere).  Bending  has  a  relatively  small  effect  on  the  circumferential 
stress. 
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Figure  12:  Stresses  at  the  junction  of  hemispherical  and  cylindrical  shells  under  external 
hydrostatic  pressure  (from  Gill,  1970.  See  also  Figure  7). 

- Inner  surface 

-  Outer  surface 
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5.6.2  Flat  Plate  End  Closure 


Although  the  hemispherical  end  closure  is  efficient  from  the  point  of  view  of 
limiting  stress  concentrations  in  a  pressure  vessel,  it  is  rather  difficult  to  fabricate,  especially 
for  large  vessels.  Rat  plate  end  closures  as  found  in  the  Type  471  submarine  design^  greatly 
simplify  construe,  The  inward  dishing  of  the  plate  under  pressure  gives  rise  to  a  large 
stress  concentration  at  the  junction  with  the  hull  cylinder.  However  the  junction  can  be 
strengthened  sufficiently  to  compensate  for  this  with  slight  weight  penalty. 

A  hemispherical  end  carries  hydrostatic  pressure  entirely  as  a  membrane  except  near  the 
junction.  In  the  case  of  a  flat  closure  the  pressure  is  carried  entirely  by  means  of  bending 
stresses.  Membrane  stresses  become  important  only  after  significant  deflection  occurs.  The 
theory  of  symmetric  bending  of  circular  plates  is  developed  by  Timoshenko  and 
Woinowsky-Kricger  (1959).  Referring  to  Figure  13  and  considering  only  axially  symmetric 
loads,  the  principal  curvatures  of  the  circular  plate  are: 

l/rn  =  -  (d2u/dr2)  =  (dp/dr)  (27) 

l/r0  =  -  (du/dr)  /  r  =  p/r  (28) 


Figure  13:  Notation  for  the  analysis  of  symmetric  bending  of  a  circular  plate. 


1  The  Type  471  end  closures  are  nol  simple  flat  plaies  because  of  adjacent  bulkheads,  tanks  and  torpedo  tubes  in  front  of 
the  forward  closure,  and  structure  supporting  the  propelior  shaft  behind  the  after  closure.  The  closures  themselves  also 
have  a  stiffening  grillage. 
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The  moments  (per  unit  length  in  the  plane  of  the  plate)  are: 

Mr  =  D(l/rn  +  v/r0)  =  D  [(d<J>/dr)  +  v<()/r]  (29) 

M0  =  D  (v/rn  +  l/r0)  =  D  [v  (d<J>/dr)  +  4>/r]  (30) 

Considering  the  torque  on  a  small  element  of  the  plate  about  a  tangential  axis  due  to  these 
moments  and  the  axially  symmetric  shear  force  Q  (per  unit  length),  the  equilibrium  condition 
is: 


Mf  +  (dMydr)r  -  M0  +  Qr  =  0 

(31) 

(d3u/dr3)  +  (d2u/dr2)  /  r  -  (du/dr)  /  r2  =  Q/D 

(32) 

In  the  case  of  hydrostatic  pressure  loading,  Q  =  pr/2  and  the  solution  of  the  differential 
equation  is: 

u  =  pr4/64D  +  0,1^74  +  02  log(r)  +  c3  (33) 


where  the  constants  c,,  c2  and  c3  are  determined  by  the  boundary  conditions.  Taking  u  =  0 
at  the  circumference  of  the  plate  (r  =  R)  and  noting  the  symmetry  requirement  (du/dr)  =  0  at 
r  =  0  we  have: 


c2  =  0 

(34) 

c3  =  -(pR4  /  64D  +  c,R2  /  4) 

(35) 

u  =  p(r4-R4)/64D  +  c1(r2-R2)/4 

(36) 

For  a  plate  loaded  with  pressure  alone,  Mr  =  0  at  r  =  R  (free  edge)  and: 

u  =  (p/64D)  (R2  -  r2)  [R2  (5  +  v)  /  (1  +  v)  -  r2]  (37) 

For  a  plate  loaded  with  a  radial  edge  bending  moment  Mr  =  MQ  (at  r  =  R)  in  addition  to  the 
pressure: 

u  =  (p/64D)  (R2  -  r2)  [R2  (5  +  v)  /  (1  +  v)  -  r2]  +  Mq(R2  -  r2)  /  2D(1  +  v)  (38) 


The  discontinuity  stresses  for  a  flat  plate  closure  of  a  cylinder  can  now  be  calculated  in  the 
same  way  as  for  a  hemispherical  cap  by  equating  deflections  and  slopes  across  the  junction. 
The  forces  required  to  maintain  continuity  arc  shown  in  Figure  14  (Cloud,  1970).  The 
response  to  pressure  and  edge  moment  have  already  been  considered.  The  shear  force  F  at 
the  edge  of  the  plate  is  required  to  balance  the  pressure  load  (F  =  pR/2)  and  is  independent  of 
the  deflection  at  the  junction.  For  plates  of  practical  proportions,  the  effect  of  the  in-planc 
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force  Qo  can  be  neglected  (the  plate  is  effectively  rigid  and  does  not  deflect  radially,  wp  =  0). 
However  F  and  Q0  must  still  be  considered  as  far  as  the  cylinder  is  concerned. 


Wa 


.M, 


i  I 


Figure  14:  Discontinuity  analysis  of  the  junction  between  a  cylindrical  shell  and  a  flat  plate 
end  closure  under  external  hydrostatic  pressure.  Forces  F  and  Qa,  and  moment  Mo  act  at  the 
junction. 


Using  equations  (4),  (17)  and  (18)  for  the  slopes  and  deflections  of  the  cylinder,  and  the 
results  above  for  the  plate,  the  continuity  conditions  at  the  junction  are: 

wp  =  0  =  wQ  +  w ,  =  2QoaR2/Eh  +  2MQa2R2/Eh  -  pR2(2  -  v)  /  2Eh  (39) 

(Note  here  that  the  term  involving  v  takes  the  end  force  F  on  the  cylinder  into  account.) 
9p  =  (du/dr)r=R  =  8o=  (dw/dz)^  (40) 

pR3/8D  (1  +  v)  +  MoR/D  (1  +  v)  =  -2Qoa2R2/Eh  -  4MQa3R2/Eh  (41) 
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(42) 


Solving  for  Mo  and  Qo  assuming  E  and  v  are  the  same  for  the  cylinder  and  the  plate: 

Mo  =  -  p  [R3a3H3  +  (2  -  v)(l  +  v)]  /4a2  (1  +  v  +  2RaH3) 

Qo  =  p  [R3a3H3  +  2  (2  -  v)(l  +  v  +  RaH3)]  /  4a(l  +  v  +  2RaH3)  (43) 

where  H  is  the  cylinderiplate  thickness  ratio  (h^)  and  a  refers  to  cylinder  dimensions. 
The  deflections  in  the  plate  can  now  be  obtained  using  equation  (38),  and  the  stresses  are 
found  as  follows  (Timoshenko  and  Woinowsky-Krieger,  1959): 

or  =  12cMr/h3  =  Ec[(d2u/dr2)  + v(du/dr)/r]  (44) 

0q  =  12cMg/h3  =  Ec[v(d2u/dr2)  +  (du/dr)  /  r]  (45) 

where  c  is  the  distance  from  the  mid-thickness  of  the  plate  and  ranges  from  -\\J2  to  \\J2. 
The  cylinder  deflections  and  stresses  due  to  the  edge  forces  Q0  and  Mo  combined  with  the 


hydrostatic  pressure  are  as  shown  below. 

w  =  2aR2  exp  (-az)  [Q0  cos(az)  +  aMo  (cos  az  -  sin  az)]  /  Eh 

-  pR2  (2  -  v)  /  2Eh  (46) 

d2w/dz2  =  4a3R2  exp  (-az)  [Qo  sin(az)  +  aMQ  (cos  az  +  sin  az)]  /  Eh  (47) 

az  =  -  pR/2h  -  Ep  [vw/R2  +  (d2w/dz2)]  /  (1  -  v2)  (48) 

a0  =  Ew/R  -  vpR/2h  -  Ep  [w/R2  +  v  (d2w/dz2)]  /  (1  -  v2)  (49) 


The  deflections  and  stresses  are  illustrated  in  Figure  15.  The  flat  plate  closure  causes  a  large 
stress  concentration  in  the  cylindrical  shell  adjacent  to  the  junction  because  of  the  bending 
induced  by  the  moment  MQ.  A  lighter  end  plate  deflects  more  under  pressure,  imposes  a 
larger  Mo  on  the  cylinder,  and  results  in  larger  stresses  in  the  cylinder.  The  stresses  in  the 
end  plate  itself  are  small  in  comparison  to  those  in  the  cylinder.  A  more  detailed  treatment  of 
this  problem  is  given  by  Watts  and  Lang  (1982).  Uddin  (1987)  considers  the  case  of  large 
deformations.  Harvey  (1985)  considers  grillage  and  ring  reinforcement  of  the  flat  plate. 

From  the  point  of  view  of  fatigue,  it  can  be  seen  that  the  flat  plate  closure  not  only 
causes  a  stress  concentration.  It  also  causes  the  stresses  near  the  outer  surface  of  the 
cylindrical  shell  to  become  tensile  near  the  junction.  In  fatigue,  only  tensile  stresses  arc 
damaging  (compressive  applied  stresses  may  also  be  damaging  in  the  presence  of  tensile 
residual  stresses). 

A  simple  flat  plate  end  closure  must  be  very  thick  compared  with  a  hemispherical  end 
closure  (or  be  reinforced).  An  alternative  would  be  to  thicken  the  cylindrical  shell  plating 
near  the  junction,  in  order  to  reduce  the  bending  stresses. 
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(a)  Circumferential  (aj  and  radial/ longitudinal  (csr  /  o2)  stresses  for  plate  thickness  500  mm. 

-  Inner  surface 

- Outer  surface 


(b)  Deflections  of  cylinder  (w)  and  plate  (u). 


-  plate  200  mm  thick 

- plate  500  mm  thick 

Figure  15:  Stresses  and  deflections  for  a  cylindrical  shell  (radius  3  m,  thickness  25  mm)  with 
a  flat  plate  end  closure  (thickness  200  mm  or  500  mm)  subjected  to  external  hydrostatic 
pressure. 
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5.6.3  Bulkhead  or  Heavy  Stiffening  Ring 

The  stress  concentration  due  to  a  bulkhead  or  heavy  stiffening  ring  may  also  be  obtained  by 
applying  the  edge  loading  formulas  in  Table  1.  Considering  an  isolated  stiffening  ring  in  an 
infinite  cylindrical  shell,  symmetry  requires  90  =  0  (Figure  10)  at  the  location  of  the  stiffener. 
This  situation  is  covered  by  Case  3  in  Table  1  (and  also  in  Table  29  of  Young,  1989). 

The  edge  force,  Q0,  of  interest  here  is  the  difference  between  the  reaction  force  of  a 
normal  stiffening  ring,  and  the  reaction  force  of  the  bulkhead  or  heavy  stiffening  ring.  The 
stress  concentration  due  to  the  bulkhead  or  heavy  stiffener  can  then  be  obtained  by 
superposition  with  the  solution  for  the  uniformly  framed  cylinder. 

The  outward  shear  force  per  unit  circumferential  length,  Qo,  applied  to  the  hull  plating  by 
a  standard  stiffening  ring  is: 

Q  =  -  Kco  0  (50) 


(The  symbol  co  refers  to  deflection  of  a  cylinder  without  a  bulkhead  and  the  symbol  w  to 
that  of  a  cylinder  with  a  bulkhead).  In  this  case  the  coordinate  origin  has  been  shifted  to  the 
location  of  the  stiffening  ring  to  simplify  the  algebra  (in  Appendices  1  to  6  the  stiffeners  are 
at  ±  L/2).  Thus  coz=0  is  the  deflection  at  the  stiffening  ring  (positive  outwards).  The 
corresponding  force  due  to  a  bulkhead  or  heavy  frame  is: 


Q1  =  -Kh^  +  w^ 


(51) 


where  the  symbols  are  defined  as  shown  in  Figure  16.  wq  is  defined  in  this  way  to  be 
compatible  with  Table  1.  Quantities,  w,  co  and  wo  are  all  positive  in  the  outward  direction,  so 
w  and  co  are  typically  negative  for  external  hydrostatic  loading. 


Cylinder 


Cylinder  +  Bulkhead 


U.F.C. 


U.F.C.  ♦  Bulkhead 


Figure  16:  Deflections  of  a  simple  cylindrical  shell  and  a  uniformly  framed  cylinder  with 
and  without  a  bulkhead. 
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The  required  edge  force  is  therefore: 

Q0  =  Q-Q 

=  (k-kh)<*Wkhwo  (52) 

From  Table  1  we  have  another  expression  for  Qo: 

Qo  =  Ehwo/R2a  (53) 


which  allows  the  following  solution  for  wq: 

wq  =  (K  -  Kh)  0)^  /  (Kh  +  Eh/R2a) 


(54) 


Again  using  Table  1  for  the  deflection  of  the  stiffener  due  to  the  bulkhead  or  heavy  stiffener, 
the  total  deflection  is: 

w  =  a)  +  R2a  (K  -  KH)  (O7^0  exp(-az)  [cos(az)  +  sin(az)]  /  (R2aK|(  +  Eh)  (55) 

where  to  is  the  deflection  of  the  uniformly  framed  cylinder  (or  of  the  simple  cylinder  in  the 
case  K  =  0). 

For  the  case  of  a  perfectly  rigid  bulkhead  (i.e.  KH  =  =  0  and  wo  =  -co)  the  hull 

plating  displacement  is: 

w  =  a)  -  (o^q  cxp(-az)  [cos(az)  +  sin(a/.)l  (56) 


At  the  other  extreme  the  bulkhead  has  the  same  stiffness  as  a  normal  stiffening  ring  (K,,  =  K) 
and  w  =  to.  Cases  intermediate  between  these  extremes  are  considered  in  Appendix  7. 

The  bulkhead  or  heavy  frame  introduces  additional  bending  stresses  compared  with  the 
uniformly  framed  cylinder  (Short  and  Bart,  1959).  For  the  worst  case  situation  of  a  perfectly 
rigid  bulkhead,  the  deflection  and  stresses  due  to  the  presence  of  the  bulkhead  arc  as  shown  in 
Figure  17.  The  stresses  are  calculated  from  the  deflection  as  in  Appendix  6  (and  combined 
with  existing  stresses  by  superposition).  The  stresses  in  the  cylinder  are  similar  to  those 
produced  by  the  flat  end  closure  (Figure  15(a)  Section  5.6.2).  Large  tensile  stresses  occur  at 
the  outer  surface  of  the  shell  plating. 


5.6.4  Other  Junctions 

Edge  bending  solutions  have  been  obtained  for  a  wide  variety  of  shells  of  revolution  by 
Hampe  (1964).  These  include  finite  as  well  as  infinite  length  cones,  cylinders  and  spherical 
segments.  Many  of  these  results  are  listed  by  Baker  ct  al.  (1972)  and  summarized  by  Young 
(1989).  The  edge  bending  solutions  allow  deflcctioas  and  stresses  to  be  calculated  for  a  wide 
variety  of  junctions  between  shells  of  revolution. 
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Figure  17:  Deflection  and  stresses  for  a  cylindrical  shell  (radius  3  m,  thickness  25  mm) 
adjacent  to  a  rigid  bulkhead. 

-  Stresses  at  inner  surface 

- Stresses  at  outer  surface 

See  also  Figures  7  and  15. 

Only  one  further  junction  will  be  considered  here;  the  junction  between  cylindrical  shells  of 
unequal  thickness.  Proceeding  as  in  Section  5.6.1,  the  deflection  at  the  jurction  (see 
Figure  18)  is. 

wi  +  wo  =  5i  +  5o  =  wo '  Pr2  (2  • v)  /  2Eh2 

=  80-pR2(2-v)/2Ehl  (57) 

80  =  -  2R2a,2  Mo  /  Eh,  -  2R2a,Qo  /  Eh,  (58) 
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W0  =  -  2R2(X22  M0  /  Eh2  +  2R2a2Q0  /  Eh2 


(59) 


O0  =  -4R2a,3  Mo  /  Ehj  -  2R2a,2Qo  /  Eh,  (60) 

=  4R2a23  Mo  /  Eh2  -  2R2022Q0  /  Eh2 


Figure  18:  Moment  Mo  and  shear  force  Qo  required  to  maintain  continuity  of  deflection  and 
slope  at  the  junction  of  two  cylindrical  shells  of  radius  R. 


Solving  the  pair  of  equations  (57)  and  (60)  gives: 
p(h,  -h2)  (2  -  v)  (ot,2^  -  cx22hl) 

M  =  -  (6 1 ) 

4(a,2h2  -  a^h,)2  -  8(a23h,  +  a,3h2)(a,h2  +  o^h,) 

Q0  =  -2Mo(a23h, +  a,3h2)/(aI2h2-a22h,) 

(62) 


Again  following  the  same  procedure  as  in  Section  5.6.1,  the  deflections  and  stresses  arc 
calculated  from  Mo  and  Qo.  Results  for  a  particular  case  are  shown  in  Figure  19.  It  should 
be  noted  that  shell  theory  does  not  take  into  account  the  stress  concentration  at  the  junction 
due  to  the  geometrical  change  of  section.  Shell  theory  assumes  uniform  distribution  of  the 
membrane  stress  and  linear  variation  of  the  bending  stress  through  the  thickness.  The  change 
of  section  leads  to  non-uniform  distribution  near  the  junction  (see  Section  5. 1 1). 
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Figure  19  :  Deflections  and  stresses  near  the  junction  of  cylindrical  shells  of  thickness 
h,  =  50  mm  and  h2  =  35  mm  (R  =  5  m,v  =  0.3)  as  calculated  from  shell  theory. 

-  Inner  surface  of  shell 

- Outer  surface  of  shell 


5.7  Non-Uniform  Framing  and  More  Complex  Geometries 

Analytical  solutions  for  submarine  hull  stresses  are  effectively  limited  to  the  simple  uniformly 
framed  cylinder  (Section  5.4)  and  to  relatively  simple  stress  concentrators  such  as  those 
considered  in  Section  5.6.  Compared  with  the  boundary  value  problem  for  the 
uniformly-framed  cylinder,  non-uniform  framing  adds  to  the  number  and  complexity  of 
boundary  conditions  to  be  satisfied.  Analysing  cylinders  of  non-uniform  thickness,  cones, 
domes,  toroids  and  combinations  of  these  forms  is  more  difficult  again  because  the  governing 
differential  equations  are  also  more  complex. 
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Kendrick  (1964)  lists  the  differential  equations  for  a  general  axisymmetric  shell.  Similar 
equations  from  Gill  (1970)  have  been  presented  in  Section  5.6  and  other  forms  may  be  found 
in  the  numerous  textbooks  on  the  subject  (e.g.  Gibson,  1980;  Kraus,  1967;  Timoshenko  and 
Woinowsky-Krieger,  1959;  Fliigge,  1960;  Baker  et  al.,  1972;  Novuzhilov,  1959).  While 
analytic  solutions  are  available  for  special  cases  of  interest  (e.g.  Young,  1989),  computer 
solutions  are  much  more  general. 

Axisymmetric  stress  analysis  is  much  simpler  than  the  general  three  dimensional 
problem,  and  computer  codes  have  been  written  to  capitalize  on  this  simplicity.  The 
programs  available  and  their  capabilities  have  been  briefly  reviewed  by  Bushnell  (1974).  His 
summary  diagram  is  shown  as  Figure  20.  A  comprehensive  review  of  the  field  is  Bushnell 
(1985).  In  addition  to  static  stress  analysis,  vibration  and  buckling  of  axisymmetric  shells  are 
encompassed  by  the  programs.  The  effects  of  axisymmetric  shape  imperfections  and  residual 
stresses  can  also  be  investigated. 

For  dealing  with  parts  of  a  submarine  structure  which  are  not  axisymmetric,  programs 
are  available  for  general  elastic-plastic  shell  analysis.  Two  of  the  better  known  ones  arc 
STAGS  (Almroth  et  al.,  1980)  and  AS  AS  (Atkins,  1982).  General  purpose  finite  element 
codes  such  as  NASTRAN  and  ABAQUS  have  also  been  applied  to  solving  design  problems. 
Brebbia  (1985)  has  surveyed  a  large  number  of  finite  element  and  boundary  element 
packages.  Their  capabilities  are  summarized  in  Figure  2 1 .  These  programs  require 
considerable  preparation  and  computing  time.  For  this  reason  they  are  only  applied  to  the 
analysis  of  structural  details  after  the  overall  structure  has  been  considered  using  an 
axisymmetric  model. 


5.8  Hull  Attachments  and  Penetrations 

We  move  now  to  analysis  of  some  non-axisymmetric  features  of  a  submarine  hull  structure. 
Local  attachments  such  as  the  fin,  and  penetrations  such  as  hatches  arc  the  most  important  of 
these.  Such  attachments  and  penetrations  usually  include  local  reinforcement  of  the  hull 
cylinder  in  the  form  of  thickened  hull  plating.  A  useful  review  of  these  aspects  is  given  by 
Kitching  (1970). 

Loads  due  to  small  attachments  and  penetrations  are  called  local  loads  because 
significant  effects  on  the  main  cylindrical  shell  are  usually  confined  to  areas  close  to  the 
attachment.  Some  typical  types  of  "branches"  are  shown  in  Figure  22.  Stress  analysis  of 
such  branches  has  been  confined  to  relatively  simple  geometries  in  which  the  boundary 
between  shell  and  attachment  has  a  relatively  simple  mathematical  description.  For 
cylindrical  shells  the  boundaries  considered  are  usually  meridional  or  circumferential  lines. 
For  spherical  shells  the  usual  boundaries  arc  circles.  Analyses  for  more  complex  shapes  are 
usually  approximate  and  arc  often  based  on  experimental  observations  utilizing  the  results  for 
the  basic  cylinder  and  sphere  intersections  (or  involve  numerical  analysis  methods). 

Under  pressure  loading,  the  change  in  shell  geometry  due  to  removing  material  for  an 
opening  or  due  to  welding  on  additional  reinforcing  material  leads  to  stress  concentration.  To 
minimize  this,  shape  changes  are  usually  gradual  (e.g.  rectangular  pads  on  a  cylinder  have 
well-rounded  comers).  Bending  stresses  can  be  expected  if  the  attachment  is  unsymmctrical 
with  respect  to  the  mid-plane  of  the  shell  plate.  In  the  following  sections,  different  types  of 
hull  penetrations  and  attachments  are  considered. 
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Figure  20:  Survey  of  computer  codes  for  stress  analysis  of  shells  and  solids  of  revolution  (Bushnell,  1974) 
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Figure  21.  Summary  of  the  capabilities  of  finite  element  and  boundary  element  computer  codes  for  stress  analysis  (Brebbia,  1985 ). 


-s- 


Figure  22:  Typical  types  of  shell  branches. 

5.8.1  Holes  in  Plates  and  Cylindrical  Shells 

Analytical  solutions  are  available  for  the  stress  concentrations  around  holes  in  flat  plates 
under  tension  or  bending,  and  also  for  simple  cylindrical  shells  under  external  pressure 
loading.  Peterson  (1974)  lists  a  variety  of  analytical  and  experimental  results  in  graphical 
form.  Consider  first  the  case  of  an  elliptical  hole  in  a  thin,  flat  plate  of  infinite  extent.  The 
classic  solution  obtained  by  Inglis  is  described  in  some  detail  by  Durelli  (1981). 

The  uniform  stresses  oz  and  oy  are  applied  to  an  infinite  sheet  containing  an  elliptical 
hole  with  semi-axes  a  and  b  as  shown  in  Figure  23(a).  The  stresses  at  the  edge  of  the  hole 
obtained  by  Inglis  are 

=  t  =  0 

(a  +  a  )  sinh  (2a  )  +  (a  -  a  )[cos(24>)  -  exp(2a  )  cos  2(<|>  -  P)] 

=  - 1 - 2 - 1 — l - - -  (63) 

cosh  (2ao)  -  cos  (2p) 


°P 


where  a  and  P  arc  elliptical  coordinates  such  that  a  =  ao  is  constant  at  the  edge  of  the 
hole  so  oa  represents  the  stress  perpendicular  to  the  edge  of  the  hole  and  Co  the  stress  tangent 
to  the  hole.  At  the  ends  of  the  major  axis  (for  4>  =  0)  z  =  a,  y  =  0,  P  =  0  ana  die  tangential 
stress  is 
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Op  =  ay(l  +2a/b)-az 


(64) 


At  the  ends  of  the  minor  axis  (<j>  =  0) 

Op  =  az  (1  +  2b/a)  -  ay  (65) 

For  the  uniaxial  stress  case  (oz  =  0)  the  stress  concentration  factor  at  the  end  of  the  major  axis 
is 


K(  =  1+2  VaTp 


(66) 


where  p  =  b2/a  is  the  radius  of  curvature  at  the  end  of  the  major  axis. 

In  the  case  of  cylindrical  bending  (as  opposed  to  biaxial  loading  above)  with  the 
geometry  shown  in  Figure  23(b),  Savin  (1961)  provides  the  following  solution.  The 
moments  at  the  edge  of  the  hole  are 

M  =  0 

r 

M.y  =  M  [1  +  2(1  +  v)(l  -  m)(m  -  cos  2y)/(3  +  v)(m2  -  2m  cos  2y+ 1)]  (67) 

where  r  and  y  are  polar  coordinates  and  m  =  (a-b)/(a+b).  The  maximum  moment  occurs  at 
y  =  ±  tc/2  and  the  minimum  at  y  =  0,  n.  The  stress  corresponding  to  the  moment  is 

o=12pMy/h3  (68) 

where  p  is  the  distance  from  the  midplane  of  the  plate  in  the  thickness  direction. 

The  stress  concentration  factor  for  a  circular  hole  in  a  cylindrical  shell  is  higher  than  for 
a  hole  in  a  flat  plate.  The  Lurd  solution  described  in  Savin’s  (1961)  compendium  is 
applicable  to  a  hole  in  a  capped,  simple  cylindrical  shell  subjected  to  pressure  p,  provided  the 
radius  of  the  hole  r  is  small  relative  to  the  shell  radius  R  and  thickness  h 

r2 «  Rh  (69) 


The  hole  is  subjected  to  uniform  biaxial  stress  oz  =  pR/2h  and  =  pR/h  as  described  in 
Section  3.  The  tangential  stress  at  the  edge  of  the  hole  is  ' 

Gy  =  (pR/h)  (3/2  +  cos(2y)  +  nr2  ^3(1  -  v2)  [4  +  5  cos(2y)]  /  8Rh)  (70) 

where  y  is  the  polar  angular  coordinate  (Figure  23b).  The  maximum  stress  occurs  at  y  =  0. 
The  stress  concentration  factor  is 


K,  =  5/2  +  9m2  V3(l  -  v2)  /  8Rh 


(71) 
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compared  with  5/2  for  a  circular  hole  in  a  flat  plate.  For  the  particular  case  R  =  3  m, 
h  =  25  mm,  r  =  250  mm  and  v  =  0.3,  the  SCF  is  7.4  (this  barely  satisfies  the  condition 
r2  «  Rh).  For  r  =  50  mm  the  SCF  is  2.7. 

Eringen  et  al.  (1965)  criticize  the  Lur6  solution  on  the  grounds  of  inappropriate  boundary 
conditions.  They  present  results  in  graphical  form  for  a  thin,  circular  cylindrical,  elastic  shell 
with  a  circular  hole,  subjected  to  various  loadings.  The  effect  of  pressure  acting  as  a  shear 
loading  on  the  edge  of  the  hole  is  taken  into  account. 

In  the  case  of  a  submarine  pressure  hull,  the  stress  around  hull  penetrations  is  reduced  by 
the  addition  of  stiffening  pads.  The  hull  stiffening  rings  also  reduce  the  circumferential 
stress  to  below  pR/h  (also  introduce  bending  moments)  but  are  not  expected  to  have  a  large 
effect  on  the  stress  concentration  factor. 


5.8.2  Local  Loads  Applied  to  Cylindrical  Shells 

Local  loads  have  generally  been  considered  in  the  form  of  line  loads,  area  loads,  or  loads 
applied  via  nozzles  or  other  local  attachments.  Summaries  are  given  by  Kitching  (1970)  and 
Kraus  (1967)  for  the  case  of  local  loading  of  cylindrical  shells.  Only  geometrically  simple 
cases  are  amenable  to  analytical  calculations. 

Hoff  and  co-workers  (Hoff,  1954;  Hoff  et  al.  1954)  obtained  analytical  results  for 
moments  and  radial  loads  applied  along  axial  lines  or  circumferential  arcs.  Expressions  for 
the  shell  deflections  u,  v  and  w  are  obtained  (containing  eight  arbitrary  constants)  which 
automatically  satisfy  the  Donnell  cylindrical  shell  equations.  The  expressions  have  the  form 

w  =  cos(nz/R) { exp(-otj 0)  [A,  cos(P,0)  +  A ^  sin(p,0)]  +  expf-a^jlAj  cos(P20)  +  A4  sin(p20)] 

+  exp(aj0)[A5  cos(Pj0)  +  Afi  sin(P[0)]  +  exp(a20)[A?  cos(P20)  +  Ag  sin(P20)]}  (72) 

Since  the  deflection  must  decay  with  distance  from  the  line  load,  only  the  negative 
exponentials  are  retained.  The  four  remaining  constants  Aj  to  A4  are  determined  by  the 
boundary  conditions.  The  stresses  and  stress  resultants  are  readily  calculated  once  the 
deflections  are  found.  Different  types  of  line  load  distributions  can  be  handled  using  Fourier 
series  methods.  The  approach  is  analogous  to  that  used  to  obtain  edge  loading  solutions  in 
Section  5.6. 

Kitching  also  aummarizes  the  analysis  of  loads  applied  over  rectangular  areas  of  a 
cylindrical  shell.  This  can  be  done  by  direct  integration  of  the  effects  of  line  loads,  or  can  be 
done  more  directly  using  Bijlaard’s  (1955a)  procedure  (see  also  Bijlaard,  1954,  1955b, 

1959a).  The  British  Standard  BS  1515  :  Pt  1  :  1965  (as  amended)  Appendix  A  contains 
numerous  graphs  allowing  calculation  of  stresses  in  cylindrical  and  spherical  shells  adjacent 
to  moments  and  radial  loads  applied  to  rectangular  areas,  axial  lines  and  circumferential  arcs. 
Some  other  relevant  documents  are  Findlay  and  Timmins  (1987)  which  gives  a  brief  review 
of  the  development  of  the  field,  British  Standard  BS  5500  Appendix  G,  and  Bedri  et  al. 

(1983).  Tooth  and  Motashar  (1989)  consider  the  case  of  radial  loading  via  a  rigid  rectangular 
attachment  rather  than  the  usual  uniformly  distributed  loading.  This  method  of  loading 
causes  stress  concentration  at  the  edges  of  the  attachment. 

The  third  major  class  of  local  loadings  which  have  been  studied  consists  of  loadings 
applied  over  more  circular  areas  of  the  cylinder  surface.  This  type  of  loading  is  important  in 
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the  case  of  intersecting  cylinders  (such  as  pressure  vessel  nozzles,  or  hatches  in  submarine 
hulls).  There  is  a  large  volume  of  literature.  Graphical  results  for  various  types  of  loading 
applicable  to  cylindrical  nozzles  in  cylindrical  pressure  vessels  are  presented  by 
Mershon  et  al.  (1984).  Gill  et  al.  (1970)  give  a  brief  overview  of  the  theory  for  branches  in 
spherical  and  cylindrical  pressure  vessels.  Teixeira  et  al.  (1981)  present  a  simplified  method 
for  calculating  stresses  due  to  loadings  on  branches  in  cylindrical  vessels.  Decock  (1973) 
and  Spence  and  Carlson  (1967)  consider  fatigue  assessment.  Leckie  and  Penny  (1963)  and 
Bijlaard  (1959b)  consider  branches  in  spherical  vessels. 

Superposition  of  stresses  due  to  external  pressure  and  those  due  to  branches  is  not 
conservative  (it  is  for  internal  pressure).  This  difficulty  can  be  overcome  using  an 
approximate  method  described  by  Bijlaard  (1955a). 


5.9  Residual  Stresses 

The  majority  of  fatigue  and  toughness  problems  associated  with  high  strength  steel  submarine 
hulls  are  limited  to  welds  and  associated  heat  affected  zones  (Roberts  and  Smith,  1988).  The 
driving  force  for  fatigue  in  these  locations  is  the  distribution  of  tensile  residual/restraint 
stresses  caused  by  welding  (Kilpatrick,  1986)  and  the  stress  concentrating  effect  of  the  weld 
profile  itself.  Without  these  stresses,  the  cyclic  compressive  stresses  due  to  hydrostatic 
loading  would  be  largely  non-damaging.  Another  significant  source  of  residual  stress  is  the 
process  of  cold  rolling  the  hull  plate  and  stiffening  rings  to  the  required  curvature.  These 
sources  arc  considered  in  Sections  5.9.1  and  5.9.2  respectively. 

In  addition  to  the  fatigue  problem,  residual  stresses  are  important  from  the  point  of  view 
of  shell  buckling.  A  brief  review  is  given  by  Bushnell  (1985,  p.  192).  Kendrick  (1985) 
points  out  that  residual  stresses  can  have  a  large  effect  on  collapse  pressure  if  initial  shape 
imperfections  are  small,  but  have  a  lesser  effect  when  the  shape  imperfections  are  greater  (i.e. 
the  two  effects  are  not  additive). 

It  is  a  relatively  simple  matter  to  calculate  the  residual  stresses  due  to  cold  bending  of 
plate.  Calculation  of  welding  stresses  is  more  difficult,  but  has  been  attempted  for  simple 
geometries  using  a  variety  of  numerical  modelling  techniques  (see  for  example  Masubuchi, 
1980;  Jones  et  al.,  1977;  Bushnell,  1985).  Only  a  simple  analysis  applicable  to  the  fatigue 
problem  is  attempted  here. 


5.9.1  Residual  Stresses  due  to  Cold  Bending  of  Plates 

Consider  the  simple  case  of  an  elastic-perfectly  plastic  plate  which  is  initially  bent  to  radius 
of  curvature  Ro  and  then  springs  back  to  radius  R.  Assuming  linear  strain  distribution 
through  the  thickness,  the  strain  and  stress  when  bent  to  radius  R0  are 


CD 

II 

1 

0 

(73) 

ae=5Ep/Ro 

1  CTy 

for  p  <  a  R0  /  E 
for  p  >  ayRQ  /  E 

(74) 
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where  p  is  the  distance  from  the  mid-plane  of  the  plate,  ay  is  its  yield  stress  and  E  its 
Young’s  modulus. 

Ignoring  Bauschinger  effect  and  assuming  elastic  springback,  the  corresponding 
quantities  when  the  plate  returns  to  radius  R  are 

e  q  =  p/R  (75) 

=  \  Ep/R  for  p  <  CT  R0  /  E  (76) 

|cy-Ep(l/R0-l/R)  forp>GyRo/E 

Equation  (76)  describes  the  residual  stress  distribution.  It  has  a  characteristic  zig-zag  shape 
as  shown  in  Figure  24. 


Figure  24:  Characteristic  residual  stress  distribution  through  the  thickness  h  of  an 
elastic-plastic  plate  after  cold  bending  to  a  curved  shape. 
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The  radius  of  curvature  R0  which  must  be  applied  to  the  plate  in  order  to  achieve  final 
radius  R  can  be  found  from  the  condition  that  the  bending  moment  must  be  zero  in  the  final 
state 


M  =  0  =  2 


Ep2/R  dp 


h/2 


+  2 


ap-Ep2(l/Ro-l/R)dp 


ayR/E 


(77) 


where  h  is  the  plate  thickness.  After  a  little  algebra  we  obtain 

Ro  =  R(4-3Q2  +  Q3)/Q3  (78) 

where  Q  is  the  parameter  Eh/cyRo.  The  quantities  defining  the  zig-zag  residual  stress 


distribution  (Figure  24)  are 

OjAry  =  R(/R  =  (4  ■  3Q2  +  Q3)  /  Q3  (79) 

a2/oy  =  (4-Q2)/  2Q2  (80) 

z/h  =  2/Q  (81) 

Kendrick  (1970)  shows  values  of  these  three  quantities  plotted  against  the  parameter 

F  =  Eh/ayR  =  R0Q/R  =  (4  -  3Q2  +  Q3)  /  Q2  (82) 


This  is  reproduced  as  Figure  25.  In  the  case  of  submarine  hulls,  F  has  a  typical  value  of  3, 
and  the  circumferential  residual  stress  peaks  due  to  cold  bending  are  about  half  the  yield 
stress.  The  outer  surface  is  in  compression  and  the  inner  surface  in  tension. 

More  sophisticated  calculations  of  the  residual  stress  distribution  have  been  made  by 
incorporating  various  work  hardening  rules  in  place  of  the  simple  elastic-perfectly  plastic 
model  considered  above  (e.g.  Shama,  1970).  However,  the  elastic-plastic  model  is  quite 
good  for  submarine  hull  steels,  and  larger  error  is  likely  to  accrue  from  lack  of  precise 
knowledge  of  the  yield  strength  of  the  particular  plate  than  from  neglect  of  work  hardening 
(or  the  Bauschinger  effect). 
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Eh  /  R°v 


Figure  25:  Effect  of  the  dimensionless  parameter  Eh/Ro>  on  the  three  parameters  defining 
the  residual  stress  distribution  (Figure  24)  of  a  plate  after  cold  bending  to  a  curved  shape. 

5.9.2  Residual  Stresses  due  to  Welding 

Welding  residual  stresses  arise  primarily  from  the  shrinkage  of  the  weld  metal  and  HAZ  on 
cooling  from  high  temperature.  The  resulting  stress  levels  depend  heavily  on  the  restraint 
imposed  on  this  shrinkage  by  the  remainder  of  the  structure.  For  high  strength  steels  typical 
of  those  used  to  construct  submarine  hulls,  a  temperature  change  of  200-300  K  is  sufficient 
for  the  shrinkage  stress  to  reach  the  yield  stress  if  thermal  movement  is  fully  restrained.  The 
residual  stress  has  short  range  and  long  range  components  due  to  non-uniform  plastic  strains 
near  the  weld  and  external  restraint  respectively  (Lidbury,  1984).  The  short  range  stresses  can 
be  considered  to  be  confined  to  distances  within  a  few  plate  thicknesses  of  the  weld.  Heavy 
section  weldments  are  normally  expected  to  have  residual  stress  close  to  the  yield  stress  in  the 
longitudinal  direction  throughout  the  fusion  zone,  and  for  a  significant  portion  of  the  joint 
thickness  in  the  transverse  direction.  The  residual  stress  in  the  through  thickness  direction  is 
much  lower  due  to  lack  of  restraint. 

Short  range  residual  stresses  are  difficult  to  treat  analytically.  However,  the  long  range 
restraint  stresses  are  more  tractable.  In  the  case  of  submarine  hulls,  circumferential  scam  and 
tec-butt  welds  have  received  the  most  attention. 
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Faulkner  (1977)  considers  two  main  weld  shrinkage  actions  in  the  tec  butt  welds  between 
the  stiffening  rings  and  hull  plate  illustrated  in  Figure  26: 

(i)  Longitudinal  (along  the  weld)  shrinkage  giving  rise  to  a  region  of  residual  stress  of 
tensile  yield  stress  magnitude  in  the  circumferential  direction.  This  region  comprises  a 
width  2  rjh  of  the  plating  and  depth  r)h  (or  0.75  T|h  according  to  Smith  and  Kirkwood, 
1977)  of  the  stiffener  web  adjacent  to  the  weld.  The  parameter  r|  varies  in  magnitude 
depending  on  the  welding  conditions,  and  the  yield  stress  and  thickness  of  the  members 
joined.  Faulkner  quotes  a  value  of  1 .5  for  a  ring  stiffened  structure.  Smith  and 
Kirkwood  quote  2-4  for  small  cylinders  with  frames  attached  by  continuous, 
single-pass  welds,  and  r)  «  2  for  submarine  hulls  and  multi-pass  welds.  The  tensile 
residual  stress  near  the  weld  is  balanced  by  an  approximately  uniform  compressive 
stress  opy,  (in  the  circumferential  direction)  distributed  over  the  remainder  of  the  cross 
section. 

(ii)  Transverse  (across  the  weld)  shrinkage  applies  a  longitudinal  bending  moment  to  the 
plating  and  causes  concertina-like  distortion.  The  additional  circumferential  stress 
op02  in  the  shell  is  balanced  by  stress  ofe  in  the  ring  stiffener. 

Although  the  longitudinal  shrinkage  would  cause  reduced  shell  diameter  at  the  ring  stiffener 
locations,  the  bending  moment  caused  by  the  transverse  shrinkage  acts  in  the  opposite 
direction.  The  net  effect  is  that  the  rings  stand  out  (for  internal  framing)  as  shown  in 
Figure  26.  The  shell  bends  inwards  between  the  rings  to  form  the  concertina  distortion 
observed  in  practice. 


_^|  2rjh  - 


Figure  26:  Residual  stresses  and  distortions  of  a  ring  stiffened  cylinder  due  to 
circumferential  tee  butt  welds  (from  Faulkner,  1977). 
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The  circumferential  stress  ap01  necessary  to  balance  the  yield  tension  block  due  to 
longitudinal  weld  shrinkage  is  easily  calculated  as 

ope,  =  -  2r)h  0^  /  (L  -  2r|h)  (83) 

where  L  is  the  spacing  between  ring  stiffeners  and  h  is  the  plate  thickness. 

Turning  to  the  transverse  shrinkage,  practical  experience  with  submarines  and  offshore 
structures  (Faulkner,  1977)  indicates  that  the  amplitude  of  the  concertina  distortion  is 
typically 


5p  =  h  / 10  (84) 

Assuming  a  sinusoidal  distortion  pattern 

8  =  8o  -  8p  sin  (nz/L)  (85) 


where  z  is  taken  to  be  zero  at  the  location  of  the  stiffening  ring  and  8o  is  yet  to  be  found 
(sec  Figure  26),  the  induced  circumferential  stresses  in  the  plate  and  in  the  frame  are 
respectively 


ap02  =  E8/R  (86) 

crf0  =  E8o/Rf  (87) 

At  equilibrium  the  loads  in  the  stiffener  and  plating  must  balance  each  other,  allowing  5o  to 
be  determined  as  follows 


0  =  ESoAf/Rf  +  2 
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Eh  [8o-8p  sin  (nz/L)]  /  R  dz 


So  =  28p  /  7t  ( l  +AfR/RfhL) 


(88) 


The  circumferential  stresses  in  the  plate  and  frame  due  to  both  longitudinal  and  transverse 
weld  shrinkage  arc  therefore 

ap0  =  f  ay  forO<z<r|h 

|  -2Tiho  /  (L-2r|h)  +  (Eh/10R)[2/n  (1  +  A/hL)  -  sin(nz/L)| 

for  Tpi  <  /.  <  L/2  (89) 

CTf0  =  Eh  /  5rcRf  (1  +  A/hL)  (90) 
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where  A  is  the  modified  frame  area  AfR/Rf.  Thu  result  reproduces  the  expressions  given  by 
Faulkner.  For  an  internally  stiffened  cylinder,  the  stiffener  is  in  tension  and  the  plating  is  in 
tension  near  the  frame  and  in  compression  near  mid-bay  as  a  result  of  the  across-lhe-weld 
shrinkage.  The  peak  compressive  stress  in  the  plating  is  of  the  order  of  a \J5  (as  measured  on 
a  submarine  hull  during  construction).  This  occurs  at  mid-bay. 

Faulkner’s  analysis  is  concerned  mainly  with  the  collapse  strength  of  ring  stiffened 
cylinders.  From  the  fatigue  viewpoint,  the  longitudinal  stresses  acting  across  the  weld  are 
more  significant  because  cracking  is  most  likely  to  initiate  at  weld  bead  toes.  The  bending 
moment  (and  hence  stresses)  can  be  calculated  from  the  curvature  of  the  plating. 

The  distortions  arising  from  welding  and  other  effects  such  as  misaligments  are 
important  not  only  because  of  the  accompanying  residual  stresses.  These  shape 
imperfections  also  modify  the  response  of  the  structure  to  externally  applied  loadings  such  as 
hydrostatic  pressure.  Buckling  pressure  can  be  greatly  reduced  by  the  presence  of  shape 
imperfections. 


5.10  Shape  Imperfections 

Shape  imperfections  can  be  broadly  classified  into  two  groups:  diffuse  imperfections  such  as 
out-of-roundness  and  concertina  distortion,  and  localised  imperfections  such  as  misfit  or 
misalignment  between  adjacent  welded  plates  or  hull  sub-units.  The  latter  group  are  better 
described  as  stress  concentrations  and  arc  considered  in  Section  5.1 1.  Generally  speaking, 
both  groups  of  shape  imperfection  introduce  additional  bending  stresses  which  would  not  be 
present  in  a  geometrically  perfect  cylindrical  or  spherical  shell. 

It  is  usual  to  consider  the  diffuse  shape  imperfections  as  having  two  components.  The 
variation  in  shape  in  the  circumferential  direction  is  known  as  out-of-roundness.  The 
variation  in  shape  in  the  longitudinal  direction  of  a  cylindrical  shell  is  known  as  the 
axisymmctric  shape  imperfection  (an  example  being  concertina  distortion  due  to  the  welded 
joints  between  hull  plating  and  ring  stiffeners).  Both  have  been  examined  in  some  detail 
because  of  the  strong  influence  they  can  have  on  the  collapse  pressure. 

Out-of-roundness  shape  imperfections  reduce  collapse  pressure  in  the  same  way  that 
eccentric  loading  of  a  column  reduces  its  buckling  load  (Wcnk,  1961;  Kendrick,  1985).  The 
degree  of  strength  reduction  depends  on  the  amplitude  and  wavelength  of  the  imperfection, 
and  also  on  the  level  of  residual  stresses.  Imperfections  arc  particularly  damaging  if  they 
match  one  of  the  buckling  modes  of  the  cylinder.  Tolerances  on  out-of-roundness  arc 
therefore  restricted  to  half  the  plate  thickness  or  less  (Wcnk,  1961). 

The  concertina  distortion  due  to  welding  has  already  been  considered  in  Section  5.9.2. 
The  initial  distortion,  comb.ned  with  the  axial  load  due  to  pressure  on  the  end  closures  of  a 
cylindrical  shell,  again  acts  to  introduce  additional  bending  stresses  as  for  an  imperfect  or 
eccentrically  loaded  column.  In  the  case  of  internal  stiffening  rings,  collapse  pressure  is 
reduced  compared  with  the  case  of  external  rings,  where  the  initial  distortion  is  outwards 
between  the  rings  (in  the  opposite  direction  to  the  displacement  under  pressure  loading). 

Fourier  scries  methods  arc  the  usual  approach  to  calculating  the  additional  stresses  due  to 
shape  imperfections.  The  initial  imperfection  can  be  expressed  in  the  form  (Kendrick,  1964): 


oo  oo 

w0  =  S  I  Cmn  cos(nO)  sin(mjtz/L) 

m=l  n=l 


(91) 


The  deflections  and  stresses  in  the  shell  due  to  hydrostatic  pressure  can  then  be  calculated. 
The  results  consist  of  a  term  representing  the  behaviour  of  a  perfect  shell,  and  an  additional 
term  due  to  the  imperfection.  Wilson  (1956b)  obtained  results  in  a  similar  form  considering 
just  the  axisymmetric  concertina  distortion.  Before  calculations  of  this  type  can  be  carried 
out,  it  is  necessary  to  know  the  actual  form  of  the  shape  imperfections. 


5.11  Stress  Concentrations 

Almost  all  of  the  discussion  in  previous  sections  has  been  concerned  with  stress  analysis 
within  the  framework  of  thin  shell  theory.  Some  of  the  basic  assumptions  of  this  theory  arc 
that  the  thickness  of  the  shell  is  small  relative  to  the  radius  of  curvature,  and  that  stresses  and 
strains  in  the  through  thickness  direction  are  small  enough  to  be  ignored.  While  the  theory 
works  well  in  describing  the  general  behaviour  of  a  structure,  it  does  not  allow  accurate 
calculation  of  stresses  at  small  scale  structural  details  where  the  geometry  changes  rapidly  on 
the  length  scale  given  by  the  thickness  of  the  shell.  However  it  is  at  just  such  structural 
details  that  fatigue  cracking  is  likely  to  initiate  because  of  stress  concentration  effects  (c.g. 
changes  of  plating  thickness,  welds).  These  higher  order  stress  concentration  factors  can  be 
calculated  if  not  too  high  (<  10)  using  a  general  finite  element  analysis  program  to  model  the 
structural  detail.  A  shell  theory  calculation  is  used  to  obtain  the  forces  and  moments  acting 
on  the  detail  at  a  distance  of  several  shell  thicknesses  away  from  the  detail  (Kendrick,  1970). 

The  numerical  approach  to  calculating  higher  order  stress  concentrations  is  beyond  the 
scope  of  this  review.  However,  a  large  body  of  information  on  stress  concentrations  has  been 
developed  by  calculation,  numerical  methods  and  experimental  methods.  Compilations  such 
as  Peterson  (1974)  provide  easy  access  to  much  of  this  information.  Some  relevant  cases 
considered  by  Peterson  arc: 

(i)  Weld  profiles  (pages  94,  95) 

(ii)  Change  of  thickness  of  cylindrical  pressure  vessel  wall  (page  97) 

(iii)  Holes  of  various  geometries  in  plate  under  a  biaxial  stress  slate  (Chapter  4) 

(iv)  Reinforced  holes  (Chapter  4) 

Additional  results  for  stress  concentration  factors  at  various  types  of  weld  arc  given  by 
Gurney  (1975,  1979a,  1979b). 

It  is  notoriously  difficult  to  predict  fatigue  behaviour  theoretically,  because  of  the 
difficulty  of  analysing  strcsscs/strains  at  structural  details.  (Further  difficulties  arise  in 
obtaining  reliable  load  histories,  and  in  taking  into  account  load  sequence  effects  and 
environmental  interactions).  The  experimental  approach  is  still  widely  used.  Research  on 
offshore  tubular  structures  has  generated  a  large  amount  of  information  available  in  the  open 
literature.  Results  for  submarine  hulls  arc  more  restricted.  One  of  the  difficulties  with  the 
experimental  approach  is  that  models  must  be  of  large  scale  in  order  to  reproduce  realistic 
welding  residual/rcstraint  stresses. 
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6.  Fatigue 


Most  of  the  information  available  on  the  fatigue  of  submarine  hulls  derives  from  the  UK 
fatigue  program  carried  out  at  ARE,  Dunfermline  (c.g.  Kilpatrick,  1986).  As  a  submarine 
hull  is  subjected  to  fluctuating  compressive  stresses,  fatigue  problems  would  not  initially  be 
expected.  However,  when  the  compressive  applied  stresses  arc  superimposed  on  tensile 
residual  welding  stresses,  at  least  part  of  the  cycle  is  tensile  and  therefore  damaging  from  a 
fatigue  viewpoint. 

The  ARE  program  includes  the  measurement  and  calculation  of  residual  and  applied 
stresses  in  scale  models  tested  to  failure  in  fatigue  chambers  (large  pressure  vessels  allowing 
pressure  cycles  to  be  applied  to  the  models).  The  work  appears  to  concentrate  on  the 
stiffener-hull  lee  butt  welds.  Allowable  fatigue  design  stresses  are  lowest  at  such 
axisymmetric  stress  concentrations.  Higher  stresses  are  allowed  at  point  stress 
concentrations  such  as  hull  penetrations.  Models  tested  to  date  failed  by  circumferential 
cracking  through  the  hull  plate,  initiating  at  the  toes  of  stiffener  tee  butt  welds  as  shown 
schematically  in  Figure  27.  Bending  stresses  after  welding  were  typically  80%  of  the  material 
yield  stress  at  the  weld  toes,  however  a  complicating  feature  of  this  type  of  work  is  that  the 
residual  stress  distribution  is  modified  during  the  first  few  load  cycles  (shakedown)  due  to 
plastic  yielding  in  compression. 


Figure  27:  Fatigue  crack  propagation  path  in  ARE  scale  model  fatigue  tests  (from 
Kilpatrick,  1986). 
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Fatigue  of  welded  structures  in  general  is  reviewed  in  great  detail  by  Gurney  (1979a). 
The  following  discussion  is  largely  derived  from  that  work.  Consider  a  residual  stress 
distribution  and  a  superimposed  cyclic  applied  stress  distribution.  If  the  yield  stress  is  not 
reached,  then  the  total  stress  is  simply  the  sum  of  the  two  components.  If  the  yield  stress  is 
reached  then  plastic  deformation  occurs,  and  that  part  of  the  structure  is  subsequently 
subjected  to  a  stress  cycle  which  ranges  from  the  yield  stress  downwards,  regardless  of  the 
stress  ratio  of  the  applied  stress.  In  many  cases  the  residual  stresses  alone  are  of  yield  stress 
magnitude  so  this  condition  is  quite  common. 

As  an  example,  consider  a  tensile  residual  stress  of  yield  stress  magnitude  (csr  =  cy)  with 
a  cyclic  applied  stress  ranging  from  -o2  to  +Oj.  The  stress  range  is  (a,  +  a2)  and  the  total 
stress  varies  between  cy  and  oy  -  (o,  +  o2).  The  fatigue  behaviour  is  dependent  only  on  the 
stress  range  and  is  independent  of  (applied)  mean  stress.  Compressive  applied  stresses  are 
just  as  damaging  as  tensile  stresses. 

Ignoring  residual  stress  redistribution,  and  considering  the  zero-compression  loading 
applicable  to  a  submarine  hull  (a,  in  the  discussion  above),  the  total  stress  ranges  from  or 
down  to  cr  -  o2.  If  o2  is  less  than  cr  then  the  stress  remains  tensile  throughout  the  cycle.  If 
o2  is  greater  than  or  then  the  total  stress  is  compressive  for  part  of  the  cycle,  the  damaging 
part  of  the  cycle  is  the  stress  range  or  to  zero,  and  fatigue  life  should  be  independent  of  a., 
(unless  o2  is  sufficiently  large  to  cause  yielding  in  compression).  It  is  clear  that  the  fatigue 
assessment  of  a  submarine  hull  design  must  encompass  both  tensile  residual  stresses  and  the 
(usually  compressive)  applied  stresses. 


7.  Conclusion 


Analytical  solutions  are  available  for  a  large  number  of  shell  structures  with  relatively  simple 
geometries  (usually  axisymmetric)  typical  of  those  found  in  submarines  and  other  pressure 
vessels.  They  give  a  good  picture  of  the  general  stress  distribution  due  to  geometrically 
simple  load  distributions  such  as  hydrostatic  pressure.  However,  analytic  solutions  arc  not 
adequate  for  obtaining  stresses  at  many  structural  details  where  fatigue  cracking  is  likely  to 
initiate.  For  these  details  the  analytical  methods  must  be  augmented  by  experimental 
methods  or  numerical  computer  calculations.  A  further  difficulty  in  predicting  fatigue 
behaviour  is  the  accurate  calculation  or  measurement  of  residual  stresses  built  into  the 
structuie  during  fabrication.  Future  work  should  address  these  questions. 


2  One  exception  to  this  statement  is  when  yielding  in  one  area  occurs  in  compression,  resulting  m  the  reduction  of  the 
balancing  residual  tensile  stress  in  an  adjacent  area.  Fatigue  crack  lip  plasticity  also  modifies  the  local  residual  stress 
distribution. 
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w  =  -  (pR2/Eh)  [1  +  2a(KR2  -  bEh)(A_sinh(az)  sin(az)  -  A+cosh(az)  cos(az))] 
d2w/dz2  =  -  (pR2)/Eh)  [4a3(KR2  -  bEh)(A+sinh(az)  sin(az)  +  A  cosh(az)  cos(az))] 

A+  =  [sinh(aL/2)  cos(aL/2)  +  cosh(aL/2)  sin(aL/2)]  /  H 

A  =  [sinh(cxL/2)  cos(aL/2)  -  cosh(aL/2)  sin(aL/2)]  /  H 

H  =  R2Ka  [sinh(aL)  +  sin(aL))  +  2Eh  [cosh(aL)  -  cos(cxL)] 
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Frame  Rigidity  Constant: 


K  =  E(Af  +  bh)  /  R2 


Stresses: 


ozm  =  0 

Ogm  =  E  e  gm  =  Ew/R 

ozb  =  -  [Ep  /  (1  -  v2)]  (d2w/dz2) 


(p  ranges  from  -h/2  at  the  inner  surface  to  +h/2  at  the  outer  surface) 


membrane 


bending 


Appendix  2 


Von  Sanden  and  Gunther  ( Including  End  Pressure) 


DE:  El  (d4w/dz4)  +  Ehw/R2  =  -  p(l  -  v/2) 


BCs:  w(-z)  =  w(z) 

dw/dz  =  0  )  at  z  =  L/2 

2EI  (d3w/dz3)  -  pb  =  Kw  j 


Particular  Solution  of  DE:  w  =  -  pR2  (1  -  v/2)  /  Eh 


General  Solution  of  DE: 

w  =  A  sinh(az)  sin(az)  +  B  sinh(az)  cos(az)  +  C  cosh(az)  sin(az) 
+  D  cosh(az)  cos(az)  -  pR2  /  (1  -  v/2) 


Solution  of  BVP: 

w  =  -  (pR2/Eh)  [(1  -  v/2)  +  2a(KR2  -  bEh)(A_sinh(az)  sin(az)  -  A+cosh(az)  cos(az))l 
d2w/dz2  =  -  (pR2/Eh)[4a3  (KR2  -  bEh)(A+sinh(az)  sin(az)  +  A  cosh(az)  cos(az))] 


Stresses: 

azm  =  -  pR  /  2h 
CTgm  =  Ew  /  R  +  va7m 

azb  =  -  [Ep  /  (1  -  v2)]  (d2w/dz2) 


membrane 


bending 
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Appendix  3 


Viterbo 


DE:  El  (d4w/d/4)  +  Ehw/R2  =  -  p(  1  -  v/2) 


BCs:  w(-z)  =  w(z) 

dw/dz  =  0  )  at  z  =  L/2 

2EI  (d3w/dz3)  -  pb  (1  -  v/2)  =  Kw  ( 


Particular  Solution  of  DE:  w  =  -  pR2  (1  -  v/2)  /  Eh 


General  Solution  of  DE: 

As  for  Von  Sanden  and  Gunther  (Including  end  pressure) 


Solution  of  BVP: 

w  =  -  [pR2  (1  -  v/2)/Eh][l  +  2a<KR2  -  bEh)(A  sinh(az)  sin(az)  -  A+cosh(az)  cos(az))] 
d2w/dz2  =  -  [pR2(l  -  v/2)/Eh)[4a3  (KR2  -  bEh)(A+sinh(az)  sin(az)  +  A  cosh(az)  cos(az))l 

Frame  Rigidity  Constant:  K  =  E  (Af  +  bh)  /  R2 


Stresses: 

a7m  =  -  pR  /  2h 

CTq"1  =  Ew  /  R  +  vazm 

azb  =  -  [Ep  /  (1  -  v2)]  (d2w/dz2) 

V  ■  ™,b 


Appendix  4 


Salerno  and  Pulos 


DE:  El  (d4w/dz4)  +  (pR/2)  (d2w/dz2)  +  Ehw/R2  =  -  p(l  -  v/2) 

BGs:  w(-z)  =  w(z) 

dw/dz  =  0  \  at  z  =  L/2 

2EI  (d3w/dz3)  -  pb(l  -  v/2)  =  Kw  j 

Particular  Solution  of  DE:  w  =  -  pR2  (1  -  v/2)  /  Eh 

General  Solution  of  DE: 

w  =  a  cxp(X,z)  +  b  cxp(^z)  +  c  cxp(^z)  +  d  exp(X4z)  -  pR2(l  -  v/2)  /  Eh 
with  X,  the  roots  of 

X4  +  (pR/2EI)  X2  +  (hw/R2I)  =  0 
giving  X,  =  ±  {2a2  1-  p/p*  ±  V( p/p*)2  - 1  ]) m 

where  p*  =  2hE/R3a2  =  2h2E/R2  '/3(1  -  v2) 

=  critical  pressure  for  axisymmctric  elastic  buckling  of  an  unstiffened  cylinder  under 
uniform  axial  pressure. 

As  X,  =  -X2  =  {2a2  [-  p/p*  +  x/(p/p*)2  -  1  ] } 1/2 
and  iX,  =  -  X4  =  {2a2  [-  p/p*  -  V(p/p*)2  -  1  ] } 1/2 
the  general  solution  reduces  to: 

w  =  A  sinh(X,z)  +  B  cosh(XjZ)  +  C  sinhfXj/.)  +  D  coshfX^.)  -  pR2  (1  -  v/2)  /  Eh 
d2w/dz~  =  AX,2  sinh(XjZ)  +  BX,2  cosh(X,z)  +  CX^  sinh/X^z)  +  DX32  coshfX^.) 
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Solution  of  BVP: 

A  =  C  =  0 

D  =  -  B  [A.,  sinh  (X,L/2)]  /  [X3  sinh  (X3L/2)] 

p(l  -  v/2)  (KR2  -  bEh)  X,  sinh  (XX/2) 

B  = - - - - - 

Eh  [G,X-  sinh  (X3L/2)  -  G3  Xj  sinh  (X,L/2)] 

where  G,  =  2EI  X;3  sinh  (X;L/2)  -  K  cosh  (X,L/2) 

Stresses: 

Gzm  =  -  pR  /  2h 

aQm  =  Ew  /  R  +  vazm 

ozb  =  -  [Ep  /  (1  -  v2)]  d2w/dz2 

a0b  =  va7b 
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Appendix  5 


Wilson’s  Asymptotic  Method 


EI  (dVdz'*)  +  (pR/2)  (d2w/dz2)  +  Ehw  /  R2  =  -  p(l  -  v/2) 


w(-z)  =  w(z) 
dw/dz  =  0 

2EI  (d3w/dz3)  -  pb(l  -  v/2)  =  Kw 


at  z  =  L/2 


Assumed  Form  of  Solution: 


w  =  pwj  +  p  w, 


Resulting  Pair  of  BVPs  for  w,  and  w2: 


EI  (d4w./dz4)  +  Ehw.  /  R2  =  -  (1  -  v/2) 


w,  (-z)  =  w.(z) 


dw,/dz  =  0 


at  z  =  L/2 


2E1  (d3w,/dz3)  -  b(l  -  v/2)  =  Kw, 


EI  (d4w2/dz4)  +  Ehw2/R2  = 


(R/2)  (d2w,/dz2) 


w2(-z)  =  w2(z) 


dw2/dz  =  0 

2EI  (d3w2/dz3>  =  Kw2 


at  z  =  L/2 


Solutions  of  BVPs: 


w  =  pw,  is  the  linear,  Viterbo  solution. 


w;  =  C,  sinh(oz)  sin(az)  +  C,  cosh(az)  cos(otz)  +  C3z  sinh(az)  cos(az) 
+  C4z  cosh(az)  sin(az) 


Cj,  C2,  C3  and  C4  arc  determined  from  the  boundary  conditions.  Fairly 
complicated  expressions  result,  as  shown  by  Wilson  (1956a). 
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Frame  Rigidity  Constant  (K):  As  for  Wilson’s  Fourier  Series  Method. 
Stresses:  As  for  Wilson’s  Fourier  Series  Method. 
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Appendix  6 


Wilson’s  Fourier  Series  Method 


DE: 

El  (d4w/dz4)  +  (vEl/R2  +  pR/2)  d2w/dz2  +  Ehw  /  R2  =  -  p(l  -  v/2) 

BCs: 

w(-z)  =  w(z) 

dw/dz  =  0  ) 

at  z  =  L/2 

2EI  (d3w/dz3)  -  pb(l  -  v/2)  =  Kw  j 

Assumed  Form  of  Solution: 

oo 

w(z)  =  w(o)  +  (2/k)  I  w(n)  cos  (2rcnz/L) 

n=l 

K 

w(n)  =  w(2tiz/L)  cos(2jtnz/L)  dz 
o 

=  Fourier  coefficients 

Solution  of  B  VP:  See  Wilson  (1956b). 

Stresses: 

CTzm  =  -  pR  /  2h  ) 

>  membrane 

CT„m  =  Ew  /  R  +  VCT  m  \ 

\J  Z  J 

ab  =  -  Ep  [vw/R2  +  (d2w/dz2)]  /  (1  -  v2)  \ 

z  bending 

agb  =  -  Ep  [w/R2  +  v  (d2w/dz2)]  /  (1  -  v2)  ) 

(p  ranges  from  -h/2  at  the  inner  surface  to  +h/2  at  the  outer  surface  of  the  shell). 

Frame  Rigidity  Constant: 

K  =  E(Af/Rf2  +  bh/R2)  +  3H(I/Rf4  +  bh3/12R4)  ^  E(Af/Rf2  +  bh/R2) 
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Appendix  7 


Stress  Concentration  due  to  Bulkhead  or  Heavy  Frame 


Calculation  of  the  hull  plating  deflection  using  equation  (55)  in  Section  5.6.3  requires 
knowledge  of  the  rigidity  constant  KH  for  the  heavy  frame.  The  special  cases  KH  =  K  (the 
normal  ring  stiffener  rigidity  constant)  and  Kf)  =  °°  (perfectly  rigid  bulkhead)  have  been 
discussed  in  Section  5.6.3.  Expressions  presented  here  allow  KH  to  be  calculated  lor 
intermediate  cases.  The  deflection  can  then  be  calculated  using  equation  (55)  and  the 
stresses  using  expressions  in  Appendix  6. 

The  assumption  is  made  that  the  additional  deflection  due  to  the  heavy  frame  is  the  same 
for  a  uniformly  framed  cylinder  as  for  an  infinite  cylindrical  shell.  This  is  a  reasonable 
assumption  for  typical  submarine  hull  deflections  and  frame  spacings. 

For  the  case  of  frames  only  moderately  more  rigid  than  the  normal  frames,  Wilson’s 
(1956c)  expressions  for  K  in  Appendix  6  may  be  used. 

Wilson  (1956c)  also  gives  an  expression  for  deep  frames  (depth  >  0.2  x  shell  radius): 

K„  =  Eb  (tf  +  h)  /  R2  +  ER,2  /  R  [(1  +  v)p  -  (1  -  v)yR,2] 

where  tf  and  additional  parameters  are  defined  in  Figure  A7-1  and: 

P  =  ?  R,2R22[Ehw±(l  -vlR^l/S 

7  =  ±R,2[Ehw*(l+v)R2Xf]/S 

\  =  Ehff  /  a2 

S  =  hwR,2  [Ehw  *  (1  +  v)  R2Xf]  -  hwR22  [Ehw  ±  (1  -  v)  R2A.f) 


The  upper  and  lower  signs  refer  to  external  and  internal  frames  respectively. 

For  a  bulkhead  idealized  as  a  fiat  circular  plate  of  thickness  t  and  radius  R,  SokolnikolT 
(1956)  shows  that  the  radial  deflection  and  stresses  due  to  uniform  radial  pressure  p  on  the 
circumference  arc: 

wr  =  -p(l  +v)(l  -2v)r/E 

CTr  =  CTe  =  -  P 

Therefore  the  rigidity  constant  is: 

Kn  =  tE  / R(1  +  v)(l  -  2v) 
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For  the  case  t  =  25  mm,  E  =  205  GPa,  R  =  5  m  and  v  =  0.3,  the  value  of  KH  is  2  GPa  and  the 
deflection  of  the  circumference  is  wr  =  0.04  mm  for  a  pressure  of  3  MPa.  The  bulkhead  is 
almost  perfectly  rigid  (the  rigidity  constant  for  a  typical  ring  stiffener  is  of  the  order  of  10% 
of  this  value). 


f 


Figure  A7-1  Notation  used  for  the  calculation  of  the  frame  rigidity  constant.  Both 
external  and  internal  frames  are  considered. 
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